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Summary 


The millimeter-wave portion of the electromagnetic spectrum is critical for 
understanding the subcloud atmospheric structure of the Jovian planets (Jupiter, 
Saturn, Uranus, and Neptune). This research utilizes a combination of laboratory 
measurements, computer modeling, and radio astronomical observation in order to 
obtain a better understanding of the millimeter-wave spectra of the Jovian planets. 
The pressure-broadened absorption from gaseous ammonia (NH 3 ) and hydrogen 
sulfide (HjS) has been measured in the laboratory under simulated conditions for 
the Jovian atmospheres. We have developed new formalisms for computing the 
absorptivity of gaseous NH S and HjS based on our laboratory measurements. We 
have developed a radiative transfer and thermochemical model to predict the abun- 
dance and distribution of absorbing constituents in the Jovian atmospheres. We 
use the model to compute the millimeter- wave emission from the Jovian planets. 
The model utilizes the results of the laboratory measurements and is also used 
to evaluate other possible candidates for millimeter-wave absorption in the Jovian 
atmospheres. Finally, we have observed Jupiter near 1.4 mm using the Caltech 
Submillimeter Observatory (CSO) in an attempt to detect gaseous hydrogen sul- 
fide. Sulfur compounds have been identified on Io, one of Jupiter’s moons, but 
they have never been detected on any of the Jovian planets. Although we were 
not able to detect hydrogen sulfide, we were able to make a good observation of 
Jupiter’s brightness temperature at this wavelength using Mars as the calibration 
standard. This research ultimately adds to the understanding of the composition 
and cloud structure of the Jovian planets and provides clues to the origin and 
evolution of the planets and solar system. 
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CHAPTER 1 
INTRODUCTION 


Section 1.1 lays the foundation for this work as it relates to the study of the Jovian 
atmospheres. The organization of the document follows in Section 1.2. 

1.1 Background and Motivation 

i 

The Jovian planets (Jupiter, Saturn, Uranus, and Neptune) are the most massive 
planetary bodies in our solar system. They are also known as the gaseous giant 
planets because a significant fraction of their mass is contained within their atmo- 
spheres. Atmospheric features such as the great red spot on Jupiter were observed 
as early as the 17th century. However, little was known about the composition 
of the Jovian atmospheres until this century. Advances in instrumentation over 
the past few decades have allowed scientists to compile much new information 
about the atmospheric structure and composition of the Jovian planets. The use 
of instruments aboard spacecraft has enhanced this new wealth of information. 

The abundances of elements observed in the atmospheres of the four Jovian 
planets provides one of the biggest clues to the origin and evolution of the planets. 
The Jovian planets have retained much of their original atmospheres, unlike then- 
inner solar system counterparts (Venus, Earth, and Mars). The most abundant 
elements in the Jovian atmospheres (and the solar system) are hydrogen and he- 
lium. Small abundances of other elements are found to exist primarily in reduced 
forms. The heavier volatiles such as CH*, NH S , and H*S appear to be enriched 
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(relative to the sun) in the atmospheres of the giant planets (see, t.g., de Pater 
et al., 1989 and Grossman, 1990). The enrichment of heavier elements favors the 
core-instability model (see, t.g., Pollack and Bodenheimer, 1989) in which the 
core of the giant planets is formed first by solid accretion. When a critical mass 
is reached, the planet begins to rapidly accrete gas from the surrounding solar 
nebula. 

Only six minor elements have been positively detected on the giant planets. 
These elements, carbon, oxygen, nitrogen, phosphorous, germanium, and arsenic, 
form a cluster in the periodic table. Sulfur remains mysteriously absent from this 
list. The formation of clouds may deplete sulfur (in the form of gaseous H 3 S) in the 
upper atmospheres of the Jovian planets, making it difficult to detect with con- 
ventional methods. Ground and space based radio observations and experiments 
provide one of best means to extract information about the presence and abun- 
dance of absorbing constituents, including HjS, below the optically thick clouds. 
One goal of this research is to use the millimeter-wave spectrum to search for 
gaseous H 2 S on Jupiter. 

Radio occultation is an example of a space-based experiment which can pro- 
vide information about the subcloud regions of the giant planets. During an oc- 
cultation, a spacecraft travels behind a planet and transmits a stable CW signal 
through the atmosphere of the planet. This signal is refracted and attenuated as 
it passes through the atmosphere of the planet. The resulting signal is received on 
the earth. Information about the planetary atmosphere can be inferred from the 
precise measurement of the signal’s frequency shift and attenuation. For example, 
the Voyager spacecraft have been used in radio occultation experiments to retrieve 
temperature-pressure profiles for all four of the Jovian planets (see, t.g., Lindal, 
tt al., 1981). These experiments typically take place at S and X Band (between 
2.3 and 8.4 GHz). Radio occultations are limited in that a single occultation can 
provide information only at one localized area of the planet. The data from only 
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one location may not be representative of the conditions elsewhere on the planet. 

The expected arrival of the Galileo spacecraft at Jupiter in 1995 and the 
Cassini mission to Saturn in the next century will provide additional clues to the 
composition of the two closest Jovian planets. These spacecraft will drop probes 
into the atmosphere of the planets. Mass spectrometers will identify gases within 
the atmosphere. Again, this type of in-situ observation is limited in that a single 
probe gives information for only one location on the planet. 

Ground-based radio astronomy has a distinct advantage over these space- 
based experiments in that the entire planet is observable with a single radio tele- 
scope or an array of telescopes. Many observations of the emission from the giant 
planets have been made with radio telescopes at wavelengths from 1 mm to several 
meters. However, the interpretation of the observations iB still in the initial stages. 
Laboratory studies of potential absorbers are needed in order to correctly inter- 
pret the measured emission from the planets. The dearth of laboratory absorption 
measurements under planetary conditions has hampered the interpretation of the 
Jovian millimeter-wave spectra in the past. Another goal of this research is to be- 
gin a program of laboratory measurements so that the available observations can 
be correctly interpreted. The millimeter- wave region of the spectrum will then be 
an important piece of the puzzle which theorists use to piece together the origin 
and evolution of the planets. 

Another difficulty in interpreting millimeter observations of the giant planets 
is the large uncertainty in the absolute flux calibration. Accurate calibration of the 
millimeter wavelength planetary observations is critical if meaningful comparisons 
are to be made between different observations and between the observations and 
radiative transfer models. Mars is the most frequently used calibrator at these 
wavelengths. However, the uncertainty in the estimated flux from Mars is reported 
to be approximately 10% (Griffin et al., 1986). Before the millimeter-wave spectra 
of the Jovian planets are fully understood, better calibration techniques will be 
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needed. 

A long standing discrepancy between modeled and observed brightness tem- 
peratures of Jupiter at millimeter wavelengths (see, e.y., de Pater and Massie, 
1985) provided the initial motivation for this work. One of the largest uncer- 
tainties in modeling the millimeter-wave emission from the giant planets is the 
absorption coefficient of gaseous ammonia (NHj). This gas is by far the strongest 
millimeter wave opacity source on Jupiter. This work began with laboratory mea- 
surements of the millimeter-wave absorption from gaseous NH3 under simulated 
Jovian conditions. The experiments were conducted at two frequency bands where 
a number of radio astronomical observations have been made. The results of the 
experiments were incorporated into a radiative transfer model which predicts the 
radio emission from the Jovian planets. The results of the radiative transfer mod- 
eling provided the motivation for a millimeter-wave observation of Jupiter in an 
attempt to detect gaseous hydrogen sulfide (HjS). A laboratory measurement of 
the pressure-broadening effects of hydrogen and helium on hydrogen sulfide was 
designed in order to correctly interpret the observation. 

This work contributes new experimental, theoretical, and observational re- 
sults which lead to a better understanding of the millimeter-wave spectra of the 
Jovian planets. Millimeter-wave instrumentation for planetary spectroscopy is still 
in an evolutionary state and will continue to provide new and improved planetary 
observations in the future. The results of the laboratory measurements presented 
in this work will help to interpret past, present, and future observations of the 
giant planets. It is hoped that this work will stimulate future research in this 
area. 


1.2 Organization 

The scope of this research may be conveniently divided into three areas: 
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1. Laboratory measurements of the millimeter-wave absorption from gaseous 
ammonia (NHj) and hydrogen sulfide (HjS) in a simulated Jovian atmo- 
sphere 

2. Radiative transfer and thermochemical modeling of the Jovian atmospheres 

3. Dual- wavelength radio astronomical observation of Jupiter at 1.4 mm 

Chapters 2, 3, and 4 discuss each of these areas, respectively. 

Chapter 2 describes the laboratory configuration, experimental methodology, 
and results of millimeter-wave absorptivity measurements of gaseous ammonia 
(NH S ) and gaseous hydrogen sulfide (HjS) under simulated Jovian conditions. 
The results are compared with various theories used to predict the millimeter- 
wave absorptivity of these gases. 

The theory and application of a radiative transfer and thermochemical model 
are described in Chapter 3. We use a forward approach to model the emission from 
each of the giant planets. In a forward approach, the parameters of the radiative 
transfer model (t.e., the abundances and distributions of absorbing constituents 
in the planetary atmospheres) are adjusted in order to obtain a good fit to the 
observed millimeter- wave spectra. The radiative transfer model utilizes the results 
of the laboratory measurements described in Chapter 2. We also develop new 
formalisms to compute the opacity from other sources in the Jovian atmospheres 
(e.g., pressure-induced absorption and water vapor absorption). 

Chapter 4 describes the approach, analysis, and results of a dual-wavelength 
radio astronomical observation of Jupiter at 1.4 mm. Relevant aspects of the 
instrumentation and calibration are described in detail. 

A summary of the major conclusions and contributions of this work is pre- 
sented in Chapter 5. In addition, this chapter provides several suggestions for 
future research. 
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CHAPTER 2 

Laboratory Measurements of Ammonia 
(NH 3 ) and Hydrogen Sulfide (H 2 S) 
Absorption Under Simulated Jovian 

Conditions 


One of the outstanding problems in the millimeter spectroscopy of planets has been 
and continues to be the lack of adequate laboratory measurements of line shapes 
and widths of gases at relevant pressures and temperatures and with appropriate 
broadening agents. This chapter describes the laboratory apparatus, procedure, 
and results of gaseous ammonia (NH S ) and hydrogen sulfide (H 2 S) absorptivity 
measurements under simulated Jovian conditions. The implications of the results 
presented in this chapter will be explored in Chapter 3. 

2.1 Propagation in a Lossy Dielectric 

The electric and magnetic fields of forward traveling waves in a lossy dielectric 
assume the form 

E{ x) = (2.1) 

H{x) = H o r ax e~ ifi \ 


and 


( 2 . 2 ) 
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respectively, where a and 0 are known as the propagation constants; a is called 
the attenuation constant or absorption coefficient, and 0 is called the phase con- 
stant. The propagation constants are related to the related to the permittivity e 
and permeability n of the medium through which the wave travels as well as the 
frequency w. In general, the permittivity of gases is complex and assumes the form 

e = e'-je". (2.3) 


The permeability n is approximately equal to the permeability of free space n 0 . 
The propagation constants for gases in general are 


a = 


and 


0 = 


respectively. Their ratio is 


a 

mat-* 


(2.4) 


(2.5) 


( 2 . 6 ) 


The loss tangent of a gaseous medium is defined as 

c «■ 
tan* = — , 

and the quality factor in a gaseous medium (Q # ) is 

i _ c 1 
g '"tan* ~V 

For a low-loss gas, the loss tangent is much less than unity. In this case, 
Equation 2.6 reduces to 

- = — . (2.9) 

0 2d K 1 


(2.7) 


( 2 . 8 ) 


R 2 * 

^ T* 


The phase constant 0 is 


(2.10) 
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so that the absorption coefficient can be written as 

(211 > 

2.2 Measurements of Ammonia (NH 3 ) Opacity 
at Ka-band and W-band 

The absorption from gaseous ammonia strongly affects the millimeter-wave spec- 
tra of the giant planets. Ammonia is by far the largest millimeter-wave opacity 
source on Jupiter and Saturn. The opacity from gaseous a mm onia must be known 
accurately before the potential effects of other absorbing constituents can be as- 
sessed. 

Steffes and Jenkins (1987) have measured the absorption from gaseous am- 
monia between 1.38 and 18.5 cm (1.6 to 22 GHz). They have shown that to within 
experimental accuracy the absorptivity of gaseous NHj is correctly expressed by 
the modified Ben-Reuven line shape as discussed by Berge and Gulkis (1976). 
However, no laboratory absorption measurements have been made under simu- 
lated Jovian conditions at frequencies above 22 GHz (wavelengths less than 1.35 
cm). Therefore, it is not known whether the use of the modified Ben-Reuven line 
shape is appropriate for computing gaseous ammonia opacity at frequencies above 
22 GHz. 

In order to test the Ben-Reuven and other line shapes, we have measured the 
opacity of gaseous ammonia (NH S ) under simulated Jovian conditions at several 
millimeter wavelengths. We conducted one set of experiments at several Ka-band 
frequencies between 32-40 GHz (7.5-9.38 mm). We conducted a second set of 
experiments at a frequency of 94 GHz or 3.2 mm (W-band). We conducted the 
experiments using mixing ratios of hydrogen, helium, and NH S which are similar to 
those fo un d on Jupiter. Jupiter’s atmosphere is approximately 90% Hj, 10% He, 
and 0.025% or 250 parts per million (ppm) NH S . In order to measure absorption, 
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we need a higher NHj mixing ratio than that found on the Jovian planets. We 
used a mixture consisting of 88.34% hydrogen (H 2 ), 9.81% helium (He), and 1.85% 
ammonia (NHj) for the Ka-band experiment. This corresponds to an NHj partial 
pressure of 28 torr within a total gas mixture pressure of 1560 torr (2 atm). The 
temperatures and pressures used in the experiments closely resemble those found 
in Jupiter’s atmosphere at altitudes which emit Ka-band radiation. We conducted 
the experiments at a pressure of 2 atm and at a temperature of 203 K. A higher 
mixing ratio of NHj is needed in order to measure W-band absorption, because 
W-band ammonia opacity is significantly less than that at Ka-band. We used a 
mixture of 85.56% H 2 , 9.37% He, and 5.07% NHj for the W-band experiment. The 
W-band experiment took place at a temperature of 210 K in order to avoid conden- 
sation. The pressures of the W-band experiments ranged from 1 to 2 atm. These 
experiments represent the first time that the opacity of gaseous ammonia has been 
measured under simulated conditions for the Jovian atmospheres at wavelengths 
less than 1 cm. 

2.2.1 Laboratory Configuration 

Figures 2.1 and 2.2 show block diagrams of the Ka-band and W-band atmospheric 
simulators. The components of the simulators may be grouped into the three sub- 
systems: an electrical subsystem, a gaseous pressure subsystem, and a temperature 
chamber (an ultra-low temperature freezer). 

The electrical subsystem of a general millimeter-wave atmospheric simulator 
is composed of an absorption cell, a millimeter-wave source, and a millimeter- 
wave receiver. For the Ka-band experiment, the electrical subsystem is a Ka- 
band Fabry-Perot resonator, a millimeter-wave, swept oscillator (Hewlett-Packard 
8690B), and a high resolution spectrum analyzer (Tektronix 7L18). The spectrum 
analyzer provides the local oscillator (LO) for an external harmonic mixer. The 
W-band electrical system is similar to the Ka-band system. The source is a power 
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Figure 2.1: Block diagram of the Ka-band atmospheric simulator 
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Figure 2.2: Block diagram of the W-band atmospheric simulator 
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supply (Micro-Now 705B) and backward-wave oscillator (BWO) tube (Micro-Now 
Model 728 RWO 110) which act together as a millimeter-wave, swept oscillator. 
The local oscillator (LO) for the harmonic mixer is a microwave source, phased 
locked to a microwave frequency counter. The mixer combines the tenth harmonic 
of the LO with the outgoing signal from the resonator. This produces an inter- 
mediate frequency (IF) of approximately 1.5 GHz. The high resolution spectrum 
analyzer displays the IF signal. 

The Ka-band Fabry-Perot resonator (shown in Figure 2.3) consists of two 
gold plated mirrors contained in a T-shaped glass pipe. The mirrors are sepa- 
rated by a distance of approximately 20 cm. The resonator is a semi-confocal 
configuration in that one mirror is flat and the other is spherical. The resonant 
frequency can be changed by turning a micrometer connected to the spherical mir- 
ror which adjusts the mirror spacing. The quality factor (Qc) of this resonator 
is appro xim ately 8000. The W-band Fabry-Perot resonator shown in Figure 2.4 
is also semi-confocal. This resonator differs from the Ka-band resonator in that 
the flat mir ror has a much smaller radius (5 cm) than the spherical mirror (11.5 
cm). The two mirrors are separated by a distance of approximately 14 cm. This 
type of configuration yields superior focusing which results in a high quality factor 
(Qc) of the resonator (approximately 25,000). The W-band resonator is contained 
in a cross shaped glass pipe with the curved mirror resting on two fixed support 
arms. The fixed arms can be adjusted in order to change the distance between 
the mirrors without disturbing the sensitive alignment. Any adjustments to the 
mirror spacing and alignment takes place before the resonator is placed in the 
temperature chamber. 

The mirrors were originally constructed of aluminum with a thin layer of gold 
sputtered on the surface of the mirrors in order to minimize resistive losses. After 
several sets of experiments, we noted a degradation in the quality factor (Qc) of 
the resonator. An inspection of the resonator revealed that ammonia had reacted 



Figure 2.3: Sketch of the Ka-band Fabry-Perot resonator 





Figure 2.4: Sketch of the W-band Fabry-Perot resonator 
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with the alu minum . This reaction caused the gold to peel off the surface of the 
mirrors. In order to correct this problem, the mirrors were refabricated using brass 
and sputtered with gold. 

The Fabry-Perot resonator operates as a band pass filter at each resonance. 
The millimeter-wave sweep oscillator is adjusted for each resonance so that it 
sweeps through the entire frequency range affected by that resonance. We measure 
the band width of each resonance with the high resolution spectrum analyzer. 
Electromagnetic energy is coupled both to and from the resonator through twin 
irises located on the flat mirror. The irises are attached to two sections of rigid 
waveguide which are sealed with rectangular pieces of mica held in place by a 
mixture of rosin and beeswax. Flexible Ka-band waveguide and rigid W-band 
waveguide connect the resonator (inside the temperature chamber) to the sweep 
oscillator and mixer (outside the temperature chamber) through a small hole in 
the temperature chamber. 

The gaseous pressure subsystem consists of the cylinders containing the var- 
ious gases used in the experiments (H2, He, NH3, and custom mixtures), two oil 
diffusion vacuum pumps, a thermocouple vacuum gauge tube (0-27 torr), a posi- 
tive pressure gauge (0-100 PSIG), and a glass pipe which contains the resonator. 
Each of the open ends of the pipe is sealed with an O-ring sandwiched between the 
glass lip and a flat brass or aluminum plate which is bolted to an inner flange. A 
network of 3/8" stainless steel tubing and valves connects the components of the 
pressure subsystem, so that each component may be isolated from the system as 
necessary (see Figures 2.1 and 2.2). When properly secured, the system is capable 
of containing two atmospheres of pressure without detectable leakage. 

Precautions have been taken to allow for the proper ventilation of the hydro- 
gen and ammonia gas so that the experiment can take place indoors. We use a 
flammable gas detector to detect any leaks during the experiment. All gases are 
released outdoors through a vent pipe where they are safely diluted by air. 
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2.2.2 Experimental Approach 

The Fabry-Perot resonator provides the sensitivity needed to measure the millimeter- 
wave aborptivity of gaseous NH S . For a Fabry-Perot resonator or interferometer, 
the quality factor of a resonance (loaded or unloaded) is expressed by 


Qc - (2-12) 

(see, t.g., Valkenburg and Derr, 1966) where A is the wavelength of the resonance, d 
is the mirror spacing, and k represents the total loss of the resonator per reflection. 
The total loss includes coupling loss, mirror loss, diffraction loss, and any loss 
due to absorbing material in the resonator. The effective path length (EPL) for 
electromagnetic energy in the resonator is the ratio of d to k or 

EPL a (2.13) 

2jt 


The addition of a lossy material in the resonator increases k which results in a 
reduced path length. The effective path length (unloaded) for both experiments 
was approximately 10 m. 

The quality factor ( Qc ) of a resonance is equal to the ratio of the energy 
stored in the resonator to the energy lost per radian. The quality factor of a 


cavity resonance in general is 



(2.14) 


(see, t.g., Terman, 1943) where fc is the center frequency of the resonance and 
BW is its full-width half-power band width. The quality factor of a resonance 


when the cavity is filled or loaded with a test gas {Qc L ) is 




(2.15) 


(Collin, 1966) where Qc v is the quality factor of the resonance when the cavity is 
evacuated, and Q t is the quality factor of the gas. Combining Equations 2.8, 2.11, 
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and 2.15 yields 

« - ( Qc\ ~ QcJ)y cm “ l > ( 2 16 ) 

where a is the absorption coefficient in cm -1 (optical depths/cm, 4.343 dB = 1 
optical depth), and X is the wavelength (in cm) of the test signal in the NH S 
mixture. 

Equation 2.16 can be simplified if the percentage change in the center fre- 
quency of the loaded {fc L ) and unloaded {fc v ) resonances is small. For our exper- 
iments, the percentage change is approximately 0.3%. Assuming that fc L = fc v 
and utilizing Equation 2.11, Equation 2.14 can be rewritten simply as 

a^(BW Ci -BW c J- (2.17) 

c 

or 

a c* 2.096 x 10 _4 (A BW) cm -1 = 90.96(ABW) dB/km, (2.18) 

where c is the speed of light, and A BW is the change in band width of the loaded 
and unloaded resonances in MHz. 

The following procedure is used to measure the unknown quantities in Equa- 
tions 2.17 and 2.18: First, we measure the unloaded band width (BWc v ) for 
each resonance while the cell is evacuated. Next, we add the hydrogen-helium- 
ammonia mixture to the cell and measure the loaded band width (BW^) for each 
resonance. The total pressure of the gas mixture is reduced by venting, and the 
measurements are repeated. Using this approach insures that the same mixture 
is used for measurements at all frequencies. Thus, even though some uncertainty 
exists in the mixing ratio and total pressure, the uncertainty for the frequency 
dependence of the millimeter-wave absorption is due only to the accuracy limits 
of the absorptivity measurements. 

The dielectric properties of non-absorbing gases such as hydrogen and helium 
can cause changes in the apparent bandwidths of resonances. The velocity (v) 
of electromagnetic energy is dependent on the dielectric constant of the medium 
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through which it travels since v = When hydrogen and helium are added 

to an evacuated resonator, a slight change occurs in the velocity or wavelength of 
the electromagnetic radiation. This slight change in wavelength can result in 
a change in coupling and a corresponding change in the quality factor or band 
width of a resonance. Because the percentage change in band width due to the 
absorption of NH S is relatively small for our system (approximately 20%), any 
changes in band width due to the dielectric effects of hydrogen and helium may 
lead to significant errors in the absorption measurement. 

The resonator, which operates as a band pass filter, is connected to a signal 
source (the millimeter-wave sweep oscillator) and to a receiver (the high resolution 
spectrum analyzer). Coupling between the resonator and the spectrum analyzer 
or sweep oscillator causes additional energy losses thereby decreasing the quality 
factor (Qc) of the resonance. The resonator was designed with minimal coupling 
in order to maximize Qc and minimize the variations in Qc that might result 
from changes in coupling that occur when gases are introduced into the resonator. 
The changes in coupling are due to the dielectric constant or permittivity of the 
test gas mixtures and are not necessarily related to the absorptivity of the gases. 
Slight imperfections in the waveguide or irises can make the apparent Qc of the 
resonator appear to vary with the abundance of lossless gases. We will refer to 
this effect as dielectric loading. 

It is necessary to repeat the absorption measurement without the absorbing 
gas present. The last step in the experimental procedure is to measure the band 
width of each resonance in a mixture consisting of 90% hydrogen (H 2 ) and 10% 
helium (He) with no ammonia present. Since the Hj/He mixture is essentially 
transparent for the pressures and wavelengths involved, no absorption is expected. 
If any apparent absorption is detected, dielectric loading (or a change in coupling 
due to the dielectric properties of the gases) is indicated. As long as the effects of 
dielectric loading are not time variable, they can be removed by measuring Q Cv 
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and BWc v with the non-ab6orbing gases present, rather than in a vacuum. 

We used pre-mixed (custom) constituent analyzed Hj/He/NHj gas mixtures 
obtained from a local gas company (Matheson) for both the Ka-band and W-band 
experiments. The gases can also be mixed in the laboratory using a thermocou- 
ple vacuum gauge. The uncertainty in the mixing ratio using the thermocouple 
vacuum gauge (±20%) is significantly higher than the uncertainty in the custom 
mixture (±2%). Before we obtained the custom mixture, we conducted one set of 
experiments at Ka-band using a mixture obtained with the thermocouple vacuum 
gauge. 

To mix the gases with the thermocouple vacuum gauge, the chamber is first 
completely evacuated, and 28 torr of gaseous ammonia is added to the system. We 
measure the pressure of the amm onia gas with the high-accuracy thermocouple 
vacuum gauge. We calibrated our vacuum gauge with a vacuum gauge from the 
Georgia Tech Microelectronics Research Center (MRC). Because the calibration 
of our vacuum gauge is based on the specific heat of air rather than ammonia, 
a conversion factor is needed to convert the gauge reading to actual ammonia 
gas pressure. We measured this conversion factor and found it to be 0.72± 0.05. 
Therefore, a reading of 20 torr on the vacuum gauge corresponds to approximately 
14 torr of actual ammonia gas pressure. 

Because the gauge is accurate only at pressures up to 20 torr, the ammonia 
gas is added in two stages. In order to obtain a total pressure of 28 torr, ammonia 
gas is initially added to the system until the thermocouple vacuum gauge reading is 
20 torr (14 torr of gaseous ammonia). We monitor the refraction due to ammonia 
gas for this pressure with the spectrum analyzer by measuring the frequency shift 
of the 39.3 GHz resonance. This frequency shift is approximately one megahertz. 
The band width change due to the absorption from self-broadened ammonia is 
negligible during this stage. Additional ammonia gas is then added to the system 
until the total frequency shift is twice that measured for an ammonia pressure of 
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14 torr. The resulting total pressure is 28 torr. 

At these low pressures, the index of refraction (relative to unity) is pro- 
portional to the ammonia gas abundance. Therefore, the ability of the system 
to accurately measure refractivity can be used to infer the relative NHj vapor 
abundance or pressure. However, it is not possible to use this approach for the 
accurate determination of absolute NHj pressure since accurate refractivity data 
for the 7.3 to 10 mm wavelength range is not yet available. In fact, by using the 
thermocouple vacuum gauge, we have measured the density normalized refractiv- 
ity of gaseous ammonia at 39.3 GHz and it found to be approximately 6.8 x 10 -17 
N-units /molecule /cm 3 . This is approximately 6 times the value at optical wave- 
lengths. The last step in the procedure is to add 1.8 atm of hydrogen (H 2 ) and 
0.2 atm of helium (He) to the system. This results in an ammonia mixing ratio of 
0.0185 at a total pressure of 2 atm (14.7 PSIG). 

2.2.3 Experimental Uncertainties 

Uncertainties in the measurement of the absorption coefficient may be classified 
into two categories: uncertainties due to instrumental error and the uncertainty 
due to noise. The uncertainties due to instrumental error are caused by the limited 
resolution and capability of the equipment used to measure pressure, temperature, 
and resonant band width. These uncertainties have been significantly reduced so 
that they are relatively small when compared with the uncertainty due to noise. 
For instance, proper calibration of the spectrum analyzer has made the uncertainty 
in the measurement of the resonant band width and center frequency in the absence 
of noise negligible. The limited ability of the low-temperature chamber to maintain 
a constant temperature results in temperature variations of only ± 2.5%. 

The largest source of uncertainty due to instrumental error in the past has 
been associated with the mixing ratio of the gas mixture. The uncertainty in the 
ammonia mixing ratio due to the thermocouple vacuum gauge is approximately 
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± 20% which results in an NH S volume mixing ratio of 0.0185± 0.0037. Even 
though measurements at all frequencies are made with the same mixing ratio 
and the frequency dependence remains intact, a large uncertainty still remains 
in the relative amplitude of the measured absorption. This uncertainty has been 
reduced by using a pre- mixed, constituent analyzed, hydrogen-helium-ammonia 
atmosphere with a mixing ratio accuracy of better than ±2% or an NHj volume 
mixing ratio of 0.0185 ± 0.00037. 

When the gaseous ammonia is introduced into the resonator and pipes, it 
adsorbs onto metallic surfaces. When the NH S gas is evacuated from the system, 
the NH 3 desorbs from the metallic surfaces leaving behind a trace amount of the 
gas. The pungent odor of gaseous NHj was indeed present when we disassembled 
the equipment (especially the stainless steel pipes) months after the experiment. 
Spilker (1990) points out that the effects and magnitudes of adsorption and des- 
orption are not fully understood and should be considered. Since our pressure 
vessel is constructed of glass and the mirrors of gold, the end plates (brass and 
aluminum) and backs of the mirrors are the only surfaces which adsorb NHj in 
the resonator. Adsorption can be a non-trivial effect in low pressure spectroscopic 
measurements. However, our measurement are made using high pressures. There- 
fore, the fractional amounts of NH3 which adsorb and desorb should be small. 

The most significant source of uncertainty in the measurement of the absorp- 
tion coe ffi cient was due to the effects of noise in the system. This electrical noise 
is displayed by the spectrum analyzer. As a result, the measurement of the band 
width of a resonance must be accompanied by an error term which is equal to the 
width of the noise on the spectrum analyzer’s display. 

In order to reduce the effects of noise, the system sensitivity must be as 
high as possible. The system sensitivity is defined as the minimum detectable 
absorptivity. It is dependent on both the quality factor (Qc) of the resonator and 
the noise present in the system. The quality factor is inversely proportional to the 
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energy lost (per cycle) in the resonator. Therefore, reducing losses in the resonator 
increases the sensitivity of the system. The loss in a Fabry-Perot type resonator 
can be attributed to three sources (Collin, 1966): resistive loss on the surfaces of 
the mirrors, coupling loss due to energy coupling out of the resonator through the 
irises, and diffraction loss around the sides of the mirrors. 

Computation of the resistive losses resulting from the gold surface of the 
mirrors showed that (in the absence of all other losses) the quality factor of the 
Ka-band resonator should have been approximately 250, (XX). However, its actual 
quality factor was approximately 10,000. The limiting factor in the performance 
of the resonator is attributed to either coupling losses or diffraction losses. In 
order to minimize the coupling losses, adjustable irises were developed so that the 
smallest possible coupling losses would occur, while still allowing sufficient signal 
coupling in and out of the resonator. However, this yielded only slightly improved 
results. 

The major limiting factor to the system sensitivity is diffraction losses around 
the edges of the mirrors. One approach used to reduce diffraction losses involves 
the precise alignment of the mirrors. We aligned the mirrors by directing the beam 
of a helium-neon laser through the input waveguide and iris and into the resonator. 
We adjusted the spherical mirror so that the reflected beam focuses precisely on 
the output iris. We found that this increased both the signal to noise ratio and the 
Qc of the resonator and therefore increased the sensitivity of the system. We used 
this tec hni que to improve the Qc of the Ka-band system. We could not use this 
technique for the W-band experiment due to curvature in the waveguide leading to 
the resonator. Fortunately, the Q c of the W-band system was already high enough 
to provide the needed sensitivity and no major adjustments were necessary. 
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2.2.4 Theoretical Characterization of Ammonia Absorp- 
tion 

The amm onia molecule (NHj) forms a tetrahedral shape as shown in Figure 2.5. 
The nitrogen atom vibrates around a stable position along an axis perpendicular to 
the plane of the hydrogen atoms. Quantum mechanics predicts that the nitrogen 
atom will tunnel through the plane of the hydrogen atoms to a stable point on the 
other side of the plane. This type of vibrational transition is known as inversion. 
Although most vibrational transitions occur in the infrared, the ammonia inver- 
sion is slowed by a peak in the potential well of the molecule which occurs as the 
nitrogen atom passes through the plane of the hydrogen atoms. This causes the 
transition to occur at microwave frequencies (Townes and Showlow, 1955). The 
inversion is coupled with rotational transitions in the molecule producing over one 
hundred individual absorption lines centered around the main inversion frequency 
of 23 GHz (1.3 cm). Under pressure, collisions with other molecules cause broad- 
ening of the individual lines. At pressures of tens of torr, the lines are broadened 
to such an extent that they overlap and form one continuous peak, obscuring the 
individual lines. 

The rotational state of a molecule is specified by two quantum numbers: J 
represents the total angular momentum which is quantized in units of the ground 
state angular momentum so that it must take on positive values; K represents the 
projection of the total angular momentum vector onto the molecular symmetry 
axis so that K < J. In general, the absorption coefficient of gaseous ammonia 
{ a NHt) is expressed as a summation over all of the rotational states: 

CKB . M = c f; £ A(J,K)F(J,K, 1 ,U.V.) cm-*, (2.19) 

J—0 K=1 

where v is frequency, 7 is the pressure-broadened line width, and u 0 is center fre- 
quency of the (J,K) transition. The summation is evaluated for J = 0, 1, 16. 
Each combination of (J,K) corresponds to a unique absorption line which results 
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Figure 2.5: Sketch of the ammonia (NH S ) molecule 


from the inversion of ammonia molecules in the ( J ,K) rotational state. An empir- 
ical correction factor, C, can be applied to bring the theory into agreement with 
laboratory data. The line intensity A can be calculated theoretically from 

A(J, K) = 1.214 ' + K)S{K)^e*W lw *- iMJ V +l) \ (2.20) 

(e.g., Berge and Gulkis, 1976) where i/ 0 (J,K) is the center frequency of the (J,K) 
transition, 7 is the pressure broadened line width, T is the temperature in kelvins, 
Pffu a is the partial pressure of NH3 gas in bar, S(K) = 3 for K a multiple of 3, 


and S(K) = 1.5 otherwise. 

The frequency dependent part of the absorption coefficient, F(J , if,'"/, 1/, !/„), 
is known as the line shape. Several different theories have been used to describe 
the line shape of collision or pressure broadened absorption lines. The Van Vleck- 
Weisskopf (1945) is 


= (£)’ • [ ( „ _ J. + h j5 + („ + J» + (,)>!' 


( 2 . 21 ) 
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where v is frequency, u 0 is the center frequency of the (J,K) transition, and 7 is the 
line width of that transition. The Van-Vleck Weisskopf line shape is known to be 
accurate at low pressures (less than 1 atm). Zhevakin and Naumov (1963) derived 
a different line shape and found that their line shape gave better results than the 
Van Vleck- Weisskopf theory when applied to atmospheric water vapor absorption 
measurements . This line shape was also derived independently by Gross (1955) 
and is sometimes referred to as the kinetic line shape. Its spectral shape is 


/ v \ 4 v ivy 

’ * “ y (i/ 0 J - */’)’ + 4 i/ j 7 j 


( 2 . 22 ) 


Ben-Reuven (1966) derived a more comprehensive line shape which was found 
to be more accurate at higher pressures. His line shape is 


F(J,if,7, S,8,u,v 0 ) 


2 / 1 \ 2 (7 ~ f)v* ± (7 + f)[fc + jf ± 7* ~ f 2 ] , 
W (j/* - (i/ a + tf) 3 - 7 J + f 2 ]* + 4 i/ j 7* 


(2.23) 


where f is a coupling element, and £ is a pressure shift term. This line shape 
reduces to the other line shapes under certain conditions. If f = 0 and 5 = 0 , 
then the Ben-Reuven line shape is equal to the Van-Vleck Weisskopf line shape. 
If 7 = f and 6 = 0 , the Ben-Reuven line shape is equal to the kinetic line shape. 

Berge and Gulkis (1976) have developed a formalism for computing the ab- 
sorption from ammonia in a hydrogen-helium atmosphere. This formalism employs 
the Ben-Reuven (1966) line shape with the parameters defined as follows: 

l(J,K) = 2.318F ft (^)’ /J +0.79P 1I .(^) ,/> +0.75P KH .(^) % (J,Jir), (2.24) 
= 1.92ft,, (^)’ /! + 0.3ft,. (^)’ /3 + 0.49ft,*. K), (2.25) 

and 

6 = 0A5P nh ,. (2.26) 


where Pn t , Ph», and iVm, are the partial pressures of hydrogen, helium, and 
ammonia in bar, and the units of 7 , and 6 are GHz. This formalism utilizes a 
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correction factor (C) in Equation 2.19 


C = 


1.0075 + 0.0304 



(2.27) 


in order to be consistent with the laboratory results of Morris and Parsons (1970), 
who measured NH 3 absorption at 9.58 GHz in a high pressure H 2 /He atmosphere 
at room temperature. 

Recently, Spilker (1990) has derived new pressure and temperature depen- 
dencies for the Ben-Reuven line shape based on high accuracy laboratory mea- 
surements under simulated Jovian conditions from 9-18 GHz. However, the ex- 
trapolation of this formalism to different pressures and frequencies can produce 


unreliable results. 

We have modified the parameters of the Ben-Reuven line shape in order to be 
compatible with the results of this work as well as the work of Morris and Parsons 
(1970), Steffes and Jenkins (1988), and Spilker (1990). This formalism employs 
the Ben-Reuven line shape as described above with the pressure-broadened line 
width and coupling element given by 



and 



respectively, where the units of 7 and f are GHz. The pressure shift term, 5, is 
the same as that of the Berge and Gulkis (1976) formalism (Equation 2.26). The 
correction term (C) used in this formalism is equal to 1. Our formalism differs 
only from the Berge and Gulkis formalism. 

Figure 2.6 shows a graph of the four theoretical formalisms which are used 
to compute ammonia absorption from 1 mm to 10 cm. The center frequencies, i/ ot 
and self-broadened line widths, 7,, of the ammonia inversion resonances used in 
all calculations are from Poynter and Kakar (1975). 


ABSORPTION dB/km 
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Ammonia absorption 



Figure 2.6: Theoretically computed ammonia absorption for 90% Hj, 10% He, 
and 0.025% NH 3 , at 2 bars, and 200 K using the Van Vleck Weisskopf lineshape 
(WW), Ben-Reuven lineshape as per Berge and Gulkis (BR-BG), Ben-Reuven 
formalism as given in this paper (BR-JS), and Zhevakin-Naumov or Gross line- 
shape (ZN-G) 
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The absorption from the pressure-broadened submillimeter rotational lines 
of NH S is significant at millimeter wavelengths. The absorption from the submil- 
limeter lines is expressed by Equation 2.19 for J = 0, 1 with 

A{J,K) = 1.826 (2.30) 

where 

B(0, 0) = 0.070(1 - e-“- 6/r ), (2.31) 

B(1,0) =0.075(1 -e" 571/r )e- M6/r , (2.32) 

and 

B{ 1, 1) = 0.053(1 - e~ 57J/T )e~ MS/T . (2.33) 

(de Pater and Massie, 1985). We use the kinetic line shape for F{J,K,~i,v,v 0 ) 
and C — 1. 


2.2.5 Experimental Results and Interpretation 

The results of the Ka-band and W-band NH S absorptivity measurements are listed 
in Table 2.1. The theoretically-derived values for the Ben-Reuven, Van Vleck- 
Weisskopf, and kinetic line shapes are also provided in tabular form. 

The Ka-band results are presented graphically in Figure 2.7. Two sets of 
measurements are shown. One set of measurements was conducted using a pre- 
mixed, constituent analyzed gas mixture (NH S mixing ratio uncertainty = ± 2% of 
its value). The gas mixture m the other set of measurements was mixed using the 
thermocouple vacuum gauge (NH S mixing ratio uncertainty = ± 20% of its value). 
Also shown in Figure 2.7 are solid lines which represent the theoretically computed 
absorption using the Van Vleck-Weisskopf line shape (WW), the Berge and Gulkis 
(1976) modified Ben-Reuven line shape (BR-BG), our modified Ben-Reuven line 
shape (BR-JS), and the Zhevakin and Naumov or Gross line shape (ZN-G). The 
Van Vleck-Weisskopf line shape overstates the opacity of NHj by nearly a factor 
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Table 2.1: Measured and theoretical values of NHj absorption 


Freq Date Press. a meas. a a a a 

(GHz) (atm) (dB/km) ZN/G WW BR-BG BR-JS 


t32.17 6/28/88 
7/7/88 
$32.89 6/28/88 
7/7/88 

f34.32 6/28/88 

7/7/88 
f35.03 6/28/88 
7/7/88 
t35.75 6/28/88 
7/7/88 
f36.46 6/28/88 
7/7/88 
|37.17 6/28/88 
7/7/88 
$37.89 6/28/88 
7/7/88 
f37.89 6/28/88 
7/7/88 
f39.32 6/28/88 
7/7/88 

$94.0 10/20/88 

10/20/88 

$94.0 10/20/88 

10/20/88 

$94.0 10/20/88 

$94.0 10/20/88 

$94.0 10/22/88 

10/22/88 


2.0 

291± 47 

2.0 

246± 45* 

2.0 

254± 53 

2.0 

210± 54* 

2.0 

200± 34 

2.0 

146± 34* 

2.0 

146± 34 

2.0 

146± 36* 

2.0 

173± 40 

2.0 

137± 45* 

2.0 

119± 40 

2.0 

146± 34* 

2.0 

155± 40 

2.0 

146± 47* 

2.0 

164± 40 

2.0 

137± 34* 

2.0 

146± 40 

2.0 

100± 34* 

2.0 

119± 34 

2.0 

109± 34* 

2.0 

115±32* 

2.0 

109±32 

1.7 

58±30* 

1.7 

36±30 

1.3 

0±30* 

1.0 

0±30* 

2.0 

91±30* 

2.0 

91±30* 


290.4 

461.6 

266.7 

435.0 

226.7 

388.6 

209.9 

368.6 

194.5 

350.0 

180.9 

333.5 

168.6 

318.4 

157.4 

304.5 

147.4 

292.0 

138.2 

280.5 

42.3 

350.6 

30.8 

253.9 

18.2 

148.9 

10.8 

88.3 

42.3 

350.6 


337.7 

321.5 

312.3 

292.0 

269.0 

244.4 

250.7 

225.2 

234.0 

208.1 

219.1 

193.3 

205.7 

180.2 

193.5 

168.5 

182.5 

158.2 

172.5 

148.8 

117.0 

120.8 

84.7 

87.4 

49.7 

51.1 

29.4 

30.3 

117.0 

120.8 


$88.34% Hj/9.81% He/1.85% NH S , T=203K 
$85.56% Hj/9.37% He/5.07% NH S , T=210K 

‘Measurements made with premixed, constituent analyzed mixture 
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of 2, while the modified Ben-Reuven line shape overstates the opacity of NH 3 by 
an average of over 40%. At these frequencies, it is not clear which line shape (if 
any) is most appropriate. 

We have further evaluated theoretical line shapes by making use of our mea- 
surements at higher frequencies where the line shapes are more distinct. Figure 2.8 
shows the results of the 94 GHz experiment as compared to four theoretical line 
shapes. At this frequency, it is clear that neither the Van-Vleck Weisskopf line 
shape nor the Gross line shape is appropriate. In Figure 2.9, we show an example 
of laboratory measurements at centimeter wavelengths by Spilker (1990) compared 
with our new formalism of the Ben-Reuven line shape. Our new line shape pro- 
vides a good fit to laboratory measurements under a wide range of temperatures, 
pressures, and frequencies. 

Theory predicts that there should be a smooth transition from the Ben- 
Reuven line shape to the Van Vleck- Weisskopf or Gross line shape at some low 
pressure. However, it is not known at what pressure this transition occurs. The 
sensitivity of the millimeter-wave measurement apparatus used in our experiments 
is not great enough to measure the absorption of a mm onia under Jovian conditions 
at pressures near or below 1 atm. In fact, we obtained no reliable Ka-band ab- 
sorption measurements at 1 atm pressure. We obtained one reliable set of W-band 
measurements at 1 atm pressure. The results are shown in Figure 2.10. Although 
we measured no significant absorption at 1 atm, the error bars give an upper limit 
for NH S opacity. This upper limit is well below the Van Vleck- Weisskopf theory 
and more compatible with the Gross line shape at 1 atm. This is most likely due 
to reduced coupling between individual lines at 1 atm pressure. 


ABSORPTION (dB/km) 
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Ammonia absorption under Jovian conditions 



FREQUENCY (GHz) 


Figure 2.7: Measured and theoretical Ka-band absorption from gaseous NH S in 
an 88.34% Hj-9.81% He-1.85% NHj mixture at 210 K and 2 atm. o: Gases mixed 
with thermocouple vacuum gauge; *: Premixed, constituent analyzed gas mixture. 

















35 


2.3 Measurement of Hydrogen Sulfide (H 2 S) Opac- 
ity at G-band 

Hydrogen sulfide (HjS) has several strong rotational lines at millimeter wave- 
lengths (168, 217, and 300 GHz). The lines are pressure broadened by H 2 and 
He in the Jovian atmospheres. We have developed a system capable of measuring 
pressure-broadening effects of Hj and He on the J' K i — Jk-iK +1 = 2o,2 — 
line at 217 GHz (1.4 mm). The laboratory configuration, procedure, and results 
are described in the following sections. 

2.3.1 Laboratory Configuration and Procedure 

Because gaseous H 2 S is extremely opaque near the center of the rotational line 
at 216.7 GHz, the extended path length provide by a resonator is not needed. 
Hydrogen sulfide absorption can be measured with a transmission cell, a millimeter 
wave source, and a detector. A block diagram of the system used to measure 
H 2 S absorption is shown in Figure 2.11. The G-Band CW signal (~218 GHz) 
is generated by doubling a W-Band (~109 GHz) klystron tube source (Varian, 

Inc. VAT 692A2). The klystron power supply (Micro-Now Model 756) provides 
1 KHz modulation by varying the voltage on the reflector of the klystron. The 
variation in frequency using this technique is less than 0.5%. Since the pressure- 
broadened line width of HjS is several GHz wide, absolute frequency stability is not 
necessary. The modulation signal incident on the frequency doubler is monitored 
with an oscilloscope. The klystron signal is sampled with a 20 dB coupler and 
downconverted to an IF of approximately 800 MHz with a harmonic mixer. A 
microwave source phase locked to a microwave frequency counter provides the LO 
for the mixer. The frequency and stability of the IF signal is monitored with a 
high resolution spectrum analyzer. The klystron frequency can be computed from 
the precise measurement of the IF and LO frequencies using the spectrum analyzer 
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Figure 2.11: Block diagram of a transmission cell for measuring H 2 S absorption 
and frequency counter. 

High gain horn antennas are used to transmit and receive the G-band signal 
which passes through a 71 cm glass cell. We aligned the antennas with a helium- 
neon laser. A G-Band square law detector converts the received millimeter-wave 
signal to a voltage which is measured with a lock in amplifier (Ithaco-Dynatrac 
391 A). A 5 cm piece of WR-5 waveguide (f e = 168 GHz) acts as a high pass filter 
to prevent any leakage of the fundamental or first harmonic (~ 109 GHz) through 
the doubler from being detected. We used a high pass filter (/ e = 300 GHz) to 
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measure the signal level of the third harmonic (~ 327 GHz). The power from the 
third harmonic was 30 dB down from the second harmonic. Thus, the detector is 
measuring power from mostly the desired second harmonic (~ 218 GHz). 

A situation analogous to dielectric loading (described in the previous section) 
occurs in a transmission cell measurement. Reflections occur at the cell boundaries 
due to the different dielectric constants of the air outside the cell, gas mixture in 
the cell, and lenses at the cell boundary. We first measured the power or voltage 
on the detector with the HjS mixture in the cell. We then measured the power 
with a gas mixture of 90% hydrogen and 10% helium in the cell. The indices 
of refraction for the two gas mixtures (with and without H 2 S, respectively) at 
STP are approximately 1.000183 and 1.000122. Reflections occurring at the cell 
boundaries should be similar for both gas mixtures. In a less rigorous check, we 
observed no difference between the signal level measured in the H 2 /He mixture and 
in a mixture of 70% H 2 /30% air which has exactly the same index of refraction 
as the H 2 S mixture. The absorption due to the hydrogen and helium mixture is 
negligible. The attenuation due to the H 2 S mixture is computed from the ratio of 
the voltages measured in the gas mixtures with and without H 2 S. This approach 
ensures that the drop in signal level is due only to absorption and not to changes 
in reflection at the cell interfaces. 

A relatively high mixing ratio of H 2 S is needed in order to measure absorption 
in the cell. The experiments take place at ambient temperature (296 K) and at a 
total pressure of 2 atm. We used a pre-mixed, constituent analyzed gas mixture 
(Matheson) in all experiments. This mixture consists of 78.79% H 2 , 9.28% He, 
and 11.93% H 2 S. 

Because H 2 S is an extremely noxious and corrosive gas, various safety pre- 
cautions are undertaken during the experiment. We used filtered gas masks and 
safety goggles when handled the gas. All equipment coming in contact with the 
H 2 S is constructed of stainless steel, glass, or plastic. The experiment takes place 
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indoors in a well ventilated area, and the gases are vented outdoors where they 
are safely diluted by air. 

2.3.2 Experimental Uncertainties 

The main source of uncertainty in this experiment is power drift in the klystron 
source. Power and frequency drifts occur as the temperature of the klystron varies. 
We found that the klystron output power exhibited a sinusoidal drift even though 
it was mounted on a large heat sink. The drift period is substantially longer 
than the time required to make an individual measurement. By obtaining several 
measurements, we can characterize the drift and minimize this uncertainty. The 
overall uncertainty in klystron power is approximately ±7%. Other instrumental 
uncertainties include uncertainty in the measurement of temperature (±1%) and 
total pressure (±7%). The uncertainty in the mixing ratio of the gas mixture is 
±2% per stated component. The total uncertainty in the measured absorption 
coefficient is the root sum square of the individual uncertainties. 


2.3.3 Theoretical Characterization of H 2 S Absorption 

In general, the opacity from a single absorption line at millimeter wavelengths is 
expressed by 

a = N S 0 (^) B/J exp[-1.439 E A*/) cm" 1 , (2.34) 

where N is the number density in molecules/cm s , T is temperature in Kelvins, T 0 
is a normalizing temperature (296 K), S 0 is a normalized line intensity, E is the 
lower state energy level, and F[v, u ot Au) is the line shape. The line parameters 
used in the computation of HjS absorption are taken from the GEISA (Gestion et 
Etude des Spectroscopiques Atmospherique) line catalog (Chedin et of., 1982 and 
Flaud et al ., 1983). 
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The Van-Vleck Weisskopf (1945) line shape used in this calculation is 

* ^ 1 (V\\ Au Au 1 

V ' Uo ' V) ~ * W l(t/ - i/ 0 ) 8 + (A*/)* + (*/ + !/.)’ + 

where A u is the pressure-broadened line width, and all frequencies are in cm -1 . 
The pressure-broadened line width of HjS in an Hj-He atmosphere is 

T n 

Au = (— ) [Afc'HjPHa + Aua,Pu« + (2.36) 

where Aur,, Ai'h,, and Ai/h,s are the hydrogen, helium and self-broadened line 
widths of H 2 S. The temperature scaling exponent, n, has not been measured for 
HjS. We assume a value of 0.67 for n, based on the values reported for the nitrogen- 
broadened line width of H 2 0 at 183 GHz (Waters, 1976). Because our measure- 
ments are conducted at room temperature, the assumed value of n does not affect 
our results. Moreover, the value of n does not significantly affect the line width 
calculation when temperatures are extrapolated to lower temperatures occurring 
in Jupiter’s atmosphere. For example, changing the value of n to 0.3 (which we es- 
timate to be a lower limit) at 200 K (the temperature near 2 bar on Jupiter) results 
in a line width difference of less than 20%. Helminger and De Lucia (1977) have 
measured the self-broadened line widths of the ~ ^K-iK+i = 2o,j — 2i ( i line 

at 217 GHz and report a value of Ai>*H a s = 9-1° MHz/Torr (6.92 GHz/bar). Willey 
et al. (1989) have measured the helium-broadened line width of the lo.i — li.o hy- 
drogen sulfide transition at 168.8 GHz and 295 K and found it to be 1.60 MHz/Torr 
(1.22 GHz/Bar). We assume the same value of Au^» for the 2 0 ,j — 2i (1 transition. 

2.3.4 Experimental Results 

The measured HjS absorption at 2 atm and 296 K is shown in Figure 2.12. The 
solid lines represent the theoretically computed absorption for several values of 
Ai'Hj. Visual inspection of Figure 2.12 suggests a value of Aur 3 approximately 
equal to 2 db 0.5 GHz/Bar (2.6 ± 0.7 MHz/Torr). 
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Figure 2.12: Measured and theoretical G-Band absorption from gaseous HjS in 
a 78.79% H a , 9.28% He, 11.93% HjS mixture at 2 atm and 296 K. Theoretical 
absorption using the Van Vleck-Weisskopf line shape for various Avb» between 1 
and 5 GHz/bar. 
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CHAPTER 3 

Modeling of the Jovian Atmospheres 


In the last chapter, we developed formalisms for computing the absorption from 
gaseous NH S and H 2 S. In this chapter, we apply the new expressions to a radiative 
transfer and thermochemical model of the giant planet atmospheres. This chapter 
begins with an explanation of a thermochemical model which is used to predict 
the distribution of volatiles (NHj, H 2 S, HjO, and CH4) in the giant planet atmo- 
spheres. Section 3.2 follows with the theory of radiative transfer as it is used to 
compute the radio emission from the giant planets. Section 3.3 describes the pa- 
rameters of the radiative transfer model in detail. In this section, we also develop 
new expressions for computing the absorption from other opacity sources. In Sec- 
tion 3.4, we adjust the parameters of the radiative transfer and thermochemical 
models for each the four Jovian planets. Using a forward modeling approach, we 
compare the modeled emission of each planet with its observed emission. 


3.1 Thermochemical Modeling 

The composition of the Jovian clouds cannot be readily determined by present 
observational techniques. However, the composition of clouds can be predicted 
with thermo chemi cal models. Weidenschilling and Lewis (1973) developed some 
of the earliest thermochemical models of the giant planets. Their models were 
based on a solar mix of elements. Recently, thermochemical models have been 
combined with radiative transfer models and used to predict the emission from 
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the giant planet atmospheres (see, t.g., Atreya and Romani, 1987 or Briggs and 
Sacket, 1989). We will use the thermochemical models in a similar approach. We 
will estimate the distribution of cloud forming constituents in the giant planet 
atmospheres and use the distributions in a radiative transfer model to predict the 
radio emission from the giant planets. We will also use the thermochemical models 
to estimate cloud bulk densities which will be used in the computation of cloud 
opacity. 

The atmosphere in the thermochemical model is composed of discrete homo- 
geneous layers. The thermochemical model begins deep in the atmosphere where 
a starting temperature and pressure have been previously established (see Sec- 
tion 3.3.1). A step dP is taken, and the average temperature and pressure are 
computed. A check is made to see if condensation has occurred. 

The saturated vapor pressure for a single constituent over a single phase is 

ln(P) = ~ + flj + flj ln(T) + 04 P “l" (3*1) 

We use values of the coefficients a,- for liquid and solid NH 3 and H 2 O from Briggs 
and Sacket (1989). The values of ai_ 3 (a 3 _ 6 = 0) for liquid and solid CH 4 are from 
de Pater and Massie (1985). Although both NH 3 and H 2 S dissolve in aqueous H 2 0 
cloud drops, the net effect on the millimeter-wave opacity, latent heat, and depres- 
sion of the freezing point is small (see, t.g., Grossman, 1990 or Briggs and Sacket, 
1989). Therefore, we will consider only a pure HjO cloud in the thermochemi- 
cal model. The coefficients a < for solid H 2 S were derived from recent laboratory 
measurements by Kraus et at. (1989). 

Gaseous NH 3 and HjS combine to form solid ammonium hydrosulfide (NRiSH) 
in the reaction 

NH S (gas) + JTjS(gas) £ NH<SH{ solid), (3.2) 

where the equilibrium constant K is related to the saturated partial pressures of 
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NHj and HjS by 


10834 

ln(ff) = ln(PNH 8 ^ffas) — 34.151 — , 


(3.3) 


where Phk» and Pjr 9 s are the partial pressures of NHj and HjS in bar. The value 
of K could be overestimated by a factor of 3 in the worst case (de Pater et al., 
1989). Therefore, the NH 4 SH cloud could form at deeper levels in the atmosphere 
(higher pressures). Since no additional laboratory data is available at this time, 
we will use Equation 3.3. 

If condensation has occurred, the incremental change in the mixing ratio of 
the condensate is computed. The mixing ratio or molar concentration (X*) for an 


ideal gas is 

X, - S, (3.4) 

where P t is the partial pressure of gas k, and P is the mean pressure. The change 
in mixing ratio, dX k , is found by differentiating Equation 3.4: 


iX„ = jdP k - y,iP. 

(3.5) 

Substituting the Clausius-Clapeyron equation, 


dP _ LP 
dT ~ RT a ’ 

(3.6) 

into Equation 3.5 yields 


dX ‘-p£> JT - L £ dR 

(3.7) 


The incremental change in the mixing ratios of NH S and HjS resulting from the 
formation of the NH 4 SH cloud is 


, Y _ ix r_ 5 sA* i[ 12 *?l < (r_ 2 < fp]. (3.8) 

dXm - - dx B , s - + Pffjs) n ^ pi 

The saturated vapor pressure for a single constituent over a single phase is 


ln(P) = ^ + oj + a 3 ln(T) + o 4 T + a B T J . (3.9) 
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We use values of the coefficients a < for liquid and solid NHj and H 2 0 from Briggs 
and Sacket (1989). The values of aj_ 2 (a s _5 = 0) for liquid and solid CH 4 are from 
de Pater and Massie (1985). Although both NH S and H 2 S dissolve in aqueous H 2 0 
cloud drops, the net effect on the millimeter-wave opacity, latent heat, and depres- 
sion of the freezing point is small (see, e.g., Grossman, 1990 or Briggs and Sacket, 
1989). Therefore, we will consider only a pure H 2 0 cloud in the thermochemi- 
cal model. The coefficients a, for solid H 2 S were derived from recent laboratory 
measurements by Kraus et ai. (1989). 

Gaseous NH S and H 2 S combine to form solid ammonium hydrosulfide (NH 4 SH) 
in the reaction 

NH 3 { gas) + J7 2 S(gas) £ NH 4 SH{boM), (3.10) 

where the equilibrium constant K is related to the saturated partial pressures of 
NHj and H 2 S by 

10834 

ln(K) = ln(/ViH 4 Pjr a s) = 34.151 — — , (3.11) 

where Pnh s and Ph 2 s are the partial pressures of NH S and H 2 S in bar. The value 
of K could be overestimated by a factor of 3 in the worst case (de Pater tt a/., 
1989). Therefore, the NH 4 SH cloud could form at deeper levels in the atmosphere 
(higher pressures). Since no additional laboratory data is available at this time, 
we will use Equation 3.8. 

The cloud density is computed according to Weidenschilling and Lewis (1973) 

D = 100 */c m* (3.12) 

dPT R ' 

where ~P is the mean pressure of the layer in bar, T is the mean temperature in 
kelvins, dP is the difference in pressure between the top and bottom of the layer, 
dXk is the change in mixing ratio of condensate k, m* is the molecular weight of 
the k th condensate, and R is the universal gas constant. The resulting cloud bulk 
densities and vertical distributions of the cloud forming constituents are integrated 
into the radiative transfer model which is described in the following section. 



3.2 Theory of Radiative Transfer 

In this section, we develop various forms of the radiative transfer equation. The 
equations are used to compute synthetic brightness spectra of the giant planets. 
We also explore the potential effect of the oblate shape of the giant planets on 
their radio emission. 


3.2.1 The Radiative Transfer Equation 

The Planck radiation law relates the brightness radiated by a black body ( B y {T )) 
to its physical temperature (T) in kelvins and frequency (i/) in Hz as follows: 

= c’[exp(£) - 1]' (313) 

where c is the speed of light (2.998 x 10 10 cm/s), h is the Planck constant (6.63 
x 10 -34 J s) and k is the Boltzmann constant (1.38 x 10 -13 J/K). The units of 
B V {T) are watts per square meter per hertz per steradian (Wm" 2 Hz _1 sr 1 ). 

The radiative transfer geometry for a planet is shown in Figure 3.1. The 
zenith angle, 0, is defined at a given point on a constant pressure surface of the 
planet. It is the angle between the observer’s line of sight to the planet and the 
local normal to the surface at that point. The quantity n is the cosine of the zenith 
angle. The radiative transfer equation gives the total brightness emitted at the 
top of the atmosphere (z = 0) along a ray in the direction /i as follows: 

£„(!» = jH «p(- jH (*•“) 

where z is the depth as measured from the top of the planet (z = 0), and a„(z) (in 
optical depths per unit length) is the total absorption coefficient or the sum of all 
of the absorbing processes at the depth z and at frequency v. The total brightness 
can also be expressed in terms of a weighting function, ll\,(z,/i): 


/ _ . \ 


dz 
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Atmospheric layers 


Figure 3.1: Sketch of the zenith angle of a planet. 

The weighting function provides a representation of which altitudes are contribut- 
ing most to the emission at a particular frequency. Another term commonly used 
in radiative transfer modeling is the vertical optical depth, r, which is defined as 

r„(z) = f a tf (z’)dz'. (3.16) 

Jo 

Substituting Equation 3.13 in Equation 3.11 yields an alternate form of the radia- 
tive transfer equation: 

B v (T,n) = j H B„(r(r)) exp(-i)^. (3.17) 


The absorption coefficient for gaseous molecules is a function of temperature, 
pressure, and frequency. The temperature, T, and pressure, P, may be related to 
depth in the atmosphere, z, by the hydrostatic equation, 


dP = p g dz , 


(3.18) 
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and the equation of state, 


P = 


pT R 


m 


(3.19) 


which yields 

d * = & p ’ (3 - 20) 
where R is the universal gas constant (8.3143 x 10 7 erg k -1 mol -1 ), m is the mean 

molecular weight of the atmosphere, and j is the acceleration due to gravity. Using 

Equation 3.17, the absorption coefficient at frequency v can be written simply as 


a function of depth z. 

The integral over r in Equation 3.14 is evaluated by modeling the atmosphere 
as a finite number of discrete layers. The average temperature, pressure, and 
altitude are specified at each layer. The integration begins at top of the atmosphere 
(P=10 mbar, t = 0, z = 0) and is terminated when r reaches a value of r 0 = 5 
optical depths: 

B,(T,n) =ES„(r(r))exp(--)— + B„(T.) exph'.h (3.21) 

where the last term is a correction factor which accounts for the brightness con- 
tribution from altitudes below the cutoff altitude, t q . The correction assumes that 
the temperature is constant below the cutoff altitude. 

Radio astronomical observations of the planets are commonly reported in 
terms of the temperature of an equivalent blackbody that would produce the ob- 
served brightness. This temperature, T B , is called the brightness temperature. It 
is found by inverting the Planck function: 



3.2.2 Disk-averaged Brightness 

The disk-averaged brightness of a planet, B y {Tp) t results from the integration of 
Equation 3.11 or 3.14 over all angles of incidence using 

B v {T d ) = 2 f 1 B v (T,n)fidfi. 

Jo 


(3.23) 
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Combining Equations 3.14 and 3.20, the disk-averaged brightness is expressed by 
B„(Jb) = 2 £ B„(T(r))exp(- T -)dr dp. (3.24) 

Spherical Approximation 

If the shape of the planet is a sphere, the integral over p may be simplified to a 
single integral (Goodman, 1969) by making use of the exponential integrals of the 
form 



(3.25) 

Substituting y = 1/p, Equation 3.21 becomes 


B.(T b ) = 2 f™ B„(t(t),h)e,(t)<It 

(3.26) 

or 


By{T D ) = -2 j^B v (T{T),n)dE*{T). 

(3.27) 


The exponential integral is evaluated numerically. 


Oblate Grid Model 


In order to test the spherical approximation, the integral over p in Equation 3.20 
can be evaluated by constructing a grid of the projection of the planet against the 
sky. The grid is made up of a finite number of rectangular elements which form the 
oblate shape of the planet. A grid of Jupiter is shown in Figure 3.2. The grid is 
more finely spaced at the edges of the planet where the limb darkening is greatest. 
The brightness temperature is computed at each point (from Equations 3.11 or 
3.14) and weighted according to the fractional area which it covers on the projec- 
tion. The disk-averaged brightness is the weighted sum of the brightness over all 
elements: 


7s (disk — averaged) 


/ x / y T B (x,y)dx dy 
Sxfy dxdy 


(3.28) 
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Figure 3.2: Grid of Jupiter. 
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Figure 3.3: Geometry of an oblate spheroid. 

The grid coordinate system is shown in Figure 3.3. The disk of the planet is 
projected onto the (x,y) plane. The line of sight from the observer to the planet is 
in the direction z. Each (z,y,z) refers to the coordinate of a point on the surface 
of the planet that is projected onto the (x,y) plane. 

The surface or radius of a Darwin-de Sitter spheroid is a function of the 
planetocentric latitude <f>' (Zharkov and Trubitsyn, 1978): 

r ($) = 1 + bo + £ 6*P„ (cos(*)) , (3.29) 

i=l 

where P 2 , is the Legendre polynomial of degree 2t, 6*,- are coefficients which are 
related to the zonal harmonic coefficients J*-, and 

* = \ - (3.30) 

The planetographic latitude, 4> is related to the planetocentric latitude <j>' by 


<f> = + orctan( ^ r ^ ). 

v r a® J 


(3.31) 
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The planetocentric longitude, A, is 

A = arctan^^j (3.32) 


and 

z = yj {r{4>) «t'n$) 5 — x 3 . (3.33) 


The 1 bar surfaces of Jupiter and Saturn are shown in Figure 3.4. For Uranus and 

Neptune, the ellipticity is less than 2%. 

The value of $ at the surface of the planet for a point (x, y ) on the projection 

can be computed using Newton’s method 


* /■ l \ r($ n )cos($„) -y 

dr(* B )cos(* n ) -r(#.)«»(* n )* 


(3.34) 


The cosine of the angle between the line of sight (z) and the local vertical or 
normal (n) is equal to n for a sphere as well as an oblate spheroid: 


H = n • z = cos(<f>) «tn(A). 


(3.35) 


For a sphere, 4> = <f>'. 

The gravitational field of the planet is expressed in terms of radial (y r ) and 
latitudinal (<^) components (see, t.g., Lindal tt a/., 1985) 


g = u r g r + g*g+ t 


(3.36) 


where u r and u i are unit vectors in the radial and latitudinal directions, respec- 
tively. The vector components are 

0rM) = - £( 2 * + 1) (“) Pu{* in <t>) + l u7 r[ 1 - (3.37) 

r r »=i r 


and 


. .. GM ^ , /R\ 2i dPn(sin<f>) 1 , dP t {sin4>) 

-J 2} 3 d<f> ’ 


(3.38) 


where G is the universal gravitational constant, M is the mass of the planet, Jj,- 
is the 2i(h zonal harmonic coefficient, R is the radius normalized to the Jj<, Pj« is 
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the Legendre polynomial of degree 2t, w is the angular velocity, r is the distance 
to the center of mass, and <f> is the planetocentric latitude. The planetocentric and 
planetographic latitudes are related to the gravitational components by 

<f> = <f>' — arctan[g+l g T ). (3.39) 

We found that the variation in gravity over the disk of the planet did not signifi- 
cantly affect the computed brightness temperature. Therefore, we use a constant 

value of g equal to the area weighted average. 

The calculation of the disk-averaged brightness temperature is computation- 
ally intensive. However, the amount of computation can be significantly reduced 
as follows: First, construct a table of the brightness temperature {T B ) as a func- 
tion of /*, from n - 0.1 to 1 (7* = 0 at /i = 0) using Equation 3.14. Then compute 
the value of n at each point on the (x,y) grid. Use a table look up with linear 
interpolation to find the brightness temperature at each point on the grid. Use a 
weighted sum to find the disk-averaged brightness temperature (Equation 3.25). 
Since the brightness as a function of n is nearly linear, the brightness temperature 
need only be computed for a few values of n (we use approximately 10). The 
symmetry of the x and y axes further decreases the amount of computation since 
disk-averaged brightness need only be computed for one quadrant of the visible 
hemisphere of the planet. 

In order to test the accuracy of the finite element grid, the total brightness 
was evaluated on a circular grid and compared with the exact solution using Equa- 
tion 3.23. The difference between the two methods was approximately 0.6% in the 
worst case. 

Finally, we computed the disk-averaged brightness with the finite element 
oblate grid and compared it with the spherical approximation. The difference 
between the two methods for both Jupiter and Saturn was less than 1% for wave- 
lengths between 1 mm and 10 cm. Since other uncertainties in the radiative trans- 
fer model are much larger than 1% (e.y., uncertainty in the absorption coefficient), 
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the use of the spherical approximation is justified. 


3.3 Parameters of the Radiative Transfer Model 
(RTM) 

The parameters of the radiative transfer model (RTM) are the temperature-pressure 
profile, absorption profiles of the absorbing constituents, and the abundance pro- 
files or vertical distributions of the absorbing constituents (abundance as a function 
of altitude or pressure). We use a forward modeling approach, in which we vary the 
parameters of the RTM in order to obtain good agreement between the modeled 
and observed spectra. 

Radiative transfer models are complex because they involve multiple param- 
eters. In most cases, a unique solution does not exist. Fortunately, some of the 
variables may be constrained. The laboratory measurements of NH3 absorption 
described in Chapter 2 are very important, because they eliminate one of the vari- 
ables in the radiative transfer model. Measurements by the Voyager spacecraft 
help to constrain other variables such as the temperature-pressure profile. In the 
following sections, we discuss the parameters of the RTM in detail. 

3.3.1 Temperature-Pressure Profile 

The computation of the temperature-pressure profile begins deep in the atmo- 
sphere of the planet where a starting temperature and pressure are specified. A 
step, dP, is taken and the corresponding change in altitude, dz is computed from 
Equation 3.17. A check is made to see if condensation has occurred. If so, the 
lapse rate is equal to the saturated adiabatic lapse rate. The saturated adiabatic 
lapse rate is found by solving the equation of energy conservation for the adiabatic 
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expansion of one mole of gas as follows: 

C p dT + ?j-dP + L k dX k = 0, (3.40) 


(see, t.g., Atreya and Romani, 1985) where T is temperature, P is pressure, R is 
the universal gas constant, dX k is the change in mixing ratio of gaseous condensate 
k, L k is the latent heat released on condensation for the k th condensate, and C p is 
the specific heat. Substituting Equations 3.7-3.8 into Equation 3.38 yields a form 
of the saturated adiabatic lapse rate: 


dT nig [l + ^ffc.LkXk + (3 , 41 ) 

iz ~ ^ [l + & (£* L ' X ‘ + ' 10845 )] ’ 

where dz is the layer thickness, m is the weighted or mean molecular weight 
of the atmosphere, g is the gravitational acceleration, and Lux is the heat of 
reaction for the formation of NH 4 SH, and C~ p is the weighted or mean specific heat 
of the atmosphere. Using the Clausius-Clapeyron equation (Equation 3.6) and 
Equation 3.1, the latent heat released per mole is expressed by 


L = + o s T + a 4 T* + 2067*, (3.42) 


where the a,- are the same coefficients as described in Section 3.1. 

If no condensation occurred, the lapse rate is equal to the dry adiabatic lapse 


rate: 

dT _ mg 
dz ~ (T p ' 

The temperature-pressure profile is computed iteratively 


(3.43) 

and constrained to meet 


the deepest level probed by the Voyager radio occultation experiments. 

A constant specific heat is used for all constituents except H,. The specific 
heats are shown in Table 3.1 (Grossman, 1990). The specific heat of hydrogen is 


explicitly dependent on the ortho-para ratio (orientation of the hydrogen atoms 
in the molecule) which is dependent on temperature. We will assume, as have 
others (see, t.g., de Pater and Massie, 1985 or Briggs and Sacket, 1989), the 
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Table 3.1: Specific heat of atmospheric constituents 


Atmospheric Constituent 

C r /R 

He (Helium) 

2.503 

CH 4 (Methane) 

4.296 

HjO (Water vapor) 

4.018 

NH S (Ammonia) 

4.459 

HjS (Hydrogen Sulfide) 

4.013 


frozen equilibrium or intermediate case for the specific heat of hydrogen (see, e.g., 
Wallace, 1980). The effect of the ortho-para ratio of hydrogen may slightly alter the 
temperature-pressure profile. However, we will consider the temperature-pressure 
profile to be a constrained variable in the radiative transfer model. 

3.3.2 Opacity 

In the last section, we constrained one of the variables of the RTM. In this section, 
we will examine the opacity from various atmospheric sources in detail. By making 
use of all of the available laboratory measurements, we will constrain this variable 
to a large extent. However, some uncertainty in the computation of opacity still 
remains. 

NHs Absorption 

The new formalism of the Ben-Reuven line shape described in Section 2.2.4 pro- 
vides the best fit to laboratory NH$ absorption measurements obtained at radio 
wavelengths. This formalism will be used to compute the absorption from gaseous 
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NH S in the radiative transfer model. 

H 2 S Absorption 

The absorption from H 2 S is computed as in Section 2.3.3. We use the line param- 
eters from the GEISA catalog (Chedin tt al., 1982) for 80 lines with frequencies 
up to 600 GHz. The computed absorption from HjS under Jovian conditions from 
1 to 600 GHz is shown in Figure 3.5. 


Pressure-Induced Absorption 


Another source of opacity in the radiative transfer model is pressure-induced ab- 
sorption from molecular hydrogen pairs H 2 -H 2 , H 2 -He, and H 2 -CH4) . Pressure- 

induced absorption is caused by the transient dipole which is induced by inter- 
molecular forces in pairs of colliding molecules. This effect is a large source of 
opacity at millimeter wavelengths for Uranus and Neptune where ammonia is de- 
pleted. 

Goodman (1969) developed a simple expression for calculating the pressure- 
induced absorption from H 2 -H 2 and H 2 -He pairs at radio wavelengths. 

4 0 x 10 -u r /273\ r8 „ _ /273 \ ,61 i 

40x10 -pj,.K(V) 1 cm • 




A 2 


(3.44) 


where Ph 2 and Ph» are the partial pressures of hydrogen and helium in atm, 
and A is the wavelength in cm. We compare this expression to the extrapolation 
of an empirical formalism developed from more recent laboratory data taken at 
infrared wavelengths (Bachet tt al ., 1983 and Dore tt al., 1983). The pressure- 
induced absorption is computed from the computationally intensive six-parameter 
empirical model given by Borysow tt al. (1985). We fit new parameters for the 
temperature and pressure dependencies in the Goodman (1969) expression to the 
Borysow tt al. (1985) model. We added a term to account for H2-CH4 absorption. 
We opt imiz ed the parameters for temperatures and pressures corresponding to 
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Figure 3.5: The computed absorption from gaseous H*S under Jovian conditions: 
XH) s 0.90, xh« = 0.10, Xj 7]5 = 3.35 x 10 P = 2 atm, T = 200 K. 
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Jupiter’s atmosphere. The new expression is 
3.557 X 10 -11 _ r_ /273 812 


s.is /273\ 2 ' 24 /273\ 884 i 

Pb*[Ph 9 {-jt) + 1.382ft, (-jr) +9.322 Pch«(-jt) ]• 

(3.45) 

Graphs of the three formalisms are shown in Figure 3.6 for 30 and 300 GHz (1 
cm and 1 mm). Our new expression deviates from the Borysow formalism by less 
than 1% for pressures less than 10 bar. At higher temperatures which correspond 
to pressures greater than 100 bar in the Jovian atmosphere, our new expression 
deviates from the Borysow formalism by approximately 10%. 


= 


A J 


HjO Absorption 

Goodman (1969) developed an expression for computing water vapor absorption 
based on the experimental data which was available at that time. This expression 
employed a Van Vleck-Weisskopf line shape for the rotational line at 22.235 GHz 
and a residual term for the remaining rotational lines which occur at higher fre- 
quencies. 

We have developed a new expression for water vapor absorption under Jo- 
vian conditions which is based on more recent laboratory data under terrestrial 
conditions. The parameters for the water vapor lines are summarized by Ulaby tt 
al. (1981) from a more detailed compilation by Waters (1976). This calculation in- 
cludes ten rotational lines with frequencies up to 448 GHz. The kinetic line shape 
is used in this calculation along with an empirical correction term derived by Gaut 
and Reifenstein (1971). Our expression is essentially the same as that described in 
Ulaby tt a/. (1981). The term P representing the pressure in the earth’s nitrogen 
and oxygen atmosphere is replaced by (0.81ft, +0.25ft e ) to reflect the broadening 
characteristics of the Jovian atmosphere. Figure 3.7 shows a comparison between 
the Goodman (1969) expression and our expression for pressures of 1 and 6 bar. 


ABSORPTION op dep/cm ABSORPTION op dep/cm 
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PRESSURE bar 



PRESSURE bar 


Figure 3.6: Computed pressure-induced absorption from Hj-Hj, Hj-He, and 
Hj-CHt pairs as a (unction of pressure Jupiter’s atmosphere at 30 and 300 GHz; 
Solid line: our expression; dot-dash line: Borysow et al. (1985) model; dashed line: 
Goodman (1969) expression. 
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Figure 3.7: Comparison between Goodman’s expression and our expression for 
HjO absorption with Xh* = Xh« = 0*1* Xff»o — 1 * 10 > A. 1 bar, 150 K, 
B: 6 bar, 300 K; Solid line: our expression, dot-dash: Goodman (1969) expression. 
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Absorption from other minor constituents 

Lellouch et al. (1984a, b) and Lellouch and Destombes (1985) have completed an 
extensive study of the detectability of minor species in the atmospheres of Jupiter 
and Saturn at millimeter wavelengths. They considered absorption from phosphine 
(PH 3 ), hydrogen cyanide (HCN), carbon monoxide (CO), and many other more 
complex molecules. The abundances of these molecules in the Jovian atmospheres 
are quite small. Emission cores arising from molecules such CO or HCN in the 
stratospheres of Jupiter or Saturn might be detectable. However, their net effect 
on the millimeter continuum is negligible. In fact, the only molecule capable of 
providing a significant amount of opacity is PH 3 . We will compute PH 3 absorp- 
tion using the parameters from the Pqynter-Pickett catalog (1981) and assuming 
the same line width as that of HjS. We will neglect the absorption from other 
molecules. 

Cloud Absorption 

Another potential source of opacity at millimeter wavelengths is scattering and 
absorption from cloud condensates. We have developed a computer program for 
computing cloud transmissivity and reflectivity using the two-stream approxima- 
tion. Our cloud model is described in detail below. In the following section, we 
use this model to compute potential cloud reflectivities and transmissivities for 
Jupiter. 

The volume extinction coefficient for cloud condensates is the sum of its 
volume absorption and scattering coefficients 

n. N. 

GCcloud = K t = K a -f- K, = ) ] Q«(^i) + ^ (3.46) 

i=l 1=1 

where Q, and Q a are the scattering and absorption cross sections, respectively, N v 
is the number of particles per unit volume, and r< is the radius of the x’** particle. 
In the Rayleigh approximation, the scattering and absorption cross sections 
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are simply 

Q. = |V|K|’ 

(3.47) 

and 

Q. = jx’lm(-K), 

(3.48) 

respectively, where 

n J + 2 ” c + 2 ’ 

(3.49) 


2xr 

X = “’ 

(3.50) 


t 

*•1 

1 

II 

<UJ 

(3.51) 


h — n — jk , 

(3.52) 


<St 

II 

<W 

(3.53) 


t' = n l - fc s , 

(3.54) 


e" = 2nk, 

(3.55) 


r is the particle radius, A is the wavelength in the same units, h is the complex 
index of refraction, and t is complex dielectric constant. The Rayleigh scattering 
approximation is acceptable for | fix |< 0*5 (Ulaby et al 1981). If | fix I-** a 
full Mie (1909) calculation should be performed. 

If the absorption cross section is much larger than the scattering cross section, 

Equation 3.44 reduces to a simple form: 


&cioud — 


18 ttM 


t 


p\ ‘■(e 1 + 2)* + e" 


si 


cm 


-l 


(3.56) 


(Battan, 1973) where p is the density of the condensation particle (in g/cm s ) , M 
is the bulk density of the cloud in the same units, A is the wavelength in cm. 

If the scattering cross section is comparable to the absorption cross section, a 
more complex model is needed. In this case, the volume absorption and scattering 


*«= [ * p[^)Qa{r)dr 

Jri 


coefficients are 


(3.57) 
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and 

K»= [ * p{r)Q,{r)dr, (3.58) 

J r i 

respectively, where r is the radius of the cloud particle, p(r) is the drop-size dis- 
tribution, and ri and r 2 are the lower and upper limits of the drop radii in the 
cloud. The phase function asymmetry, g, is 

9 = — / G p(r) Q,(r)dr, (3.59) 

K, Jr i 

where G is the fraction of power scattered in the forward direction. We assume 
that the scattering is from spherical particles, so that G is determined from the 
Mie coefficients (see e.g., Gasiewski, 1990). 

The uncertainties in modeling the cloud absorption are very large and arise 
from several sources. One uncertainty is the cloud bulk density. We estimate the 
densities with the thermochemical model described in Section 3.1. Romani (private 
communication, 1990) suggests that the actual cloud densities on the planets may 
be a factor of 5 less than those predicted with equilibrium condensation models 
based on studies of the earth’s clouds. Therefore, the cloud bulk density will be a 
variable in the computation of cloud opacity. Another uncertainty in calculating 
cloud opacity is the complex index of refraction of NHs and NH 4 SH ice. The 
complex index of refraction for NHs ice has been measured at infrared wavelengths 
by Sill et al. (1980). We use a value of 1.3 for (Sill tt al., 1980), 1.74 for 
n NffiSH (CRC, 1980), and 1.78 for ice (Ulaby, 1981). The value of n should 
not differ significantly from infrared to millimeter wavelengths. However, the value 
of k is frequency dependent. We will assume *nh, and £nh 4 sh to be variables. For 
HjO ice, 

*ir 9 o = (3.60) 

(Grossman, 1990) where a is the absorption coefficient. At A = 1 mm, k^o = 
0.00035. 

Finally, there is a large uncertainty in the cloud particle size and distribution. 
West et al. (1986) have given a summary of constraints on the particle sizes of 
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Jupiter’s clouds based on a compilation of available observations. There is evidence 
for a stable layer of particles having effective radii of 1 or smaller at pressures 
from 0.2 to 0.7 bar. A layer of larger particles (~ 3 - 100 #zm) may exist near 0.7 
bar primarily in Jupiter’s zones. This cloud is thought to be a mixture of NH S ice 
and chromophores (although there is no direct evidence). Bjoraker et al. (1986a, b) 
used a massive absorbing cloud (greybody with r ~ 3) at 2 bar in order to explain 
observations at 5 /im. This second cloud is thought to be NH4SH ice and/or HjO 
ice. There are no estimates of particle sizes for this cloud. 

We have assumed the particle-size distribution to be a modified gamma dis- 
tribution (Deirmendjian, 1969): 

p(r) = or“exp(— 6r' 1 ), (3.61) 

where a, b , a, and 7 are shape factors. We assume two of the shape factors to be 
similar to those of a cirrostratus cloud: 

a = 6, P = 0.5, ( 3 - 62 ) 


(see, e.g., Ulaby et al ., 1981). The remaining shape factors b and 7 are determined 
by specifying a mode radius (r e ) and bulk density M for the cloud. The mode 
radius, r e , is defined as the radius at which the distribution p(r) is at a maximum. 
The shape factors may be written in terms of the mode radius, 



(3.63) 


and the bulk density, M, 


M — 


pa* 
3 7 6* 


T(0) g/cm 3 , 


(3.64) 


where 


p is the density of the condensate, and T is the gamma function. 


(3.65) 
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The radiative transfer problem is simplified by assuming that radiation is 
scattered in only the forward and backward directions. The cloud is assumed to 
be made up of discrete homogeneous layers. The reflection and transmission for a 
layer of thickness dz are 


or 

i + vT — w J coth(a e dz) 


(3.66) 


and 

_ \/l — u?csch(a e dz) 

e 1 + y/l — w J coth(cr c dz) ' 
respectively (Gasiewski, 1990), where a e is the mode power attenuation coefficient 



<*e = \J*l + 2/c a *c, t , 


(3.68) 


JjJ is the back-scattering albedo 


u = 


Kth 

K a + K ti ’ 


(3.69) 


K a includes both gaseous absorption and absorption from cloud particles, and k,i 
is the back-scattering coefficient 


K th = K t 


1 -9 
2 ' 


(3.70) 


If the reflectivity and transmissivity of two adjacent layers (a and b) are specified, 
the combined reflectivity and transmissivity of both layers are 


r.* 


r* + 


1 - r a r 4 


(3.71) 


and 

= l-V.r.’ (372) 

(e.g., Paltridge and Platt, 1976) respectively, for downwelling radiation. A similar 
set of equations applies to upwelling radiation. 

Once the reflectivity and transmissivity are known, the brightness temper- 
ature at n = 1 can be roughly estimated using the following approach: Suppose 
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that there is a cloud at pressure P e and temperature T e , with reflectivity re e and 
transmissivity tr e , and that the zenith optical depth above the cloud is t x . The 
total brightness or emission at n = 1 is the sum of the emission from the atmo- 
sphere above the cloud, emission from the cloud itself, emission reflected by the 
cloud, and emission from below the cloud: 

B v {T,fi = 1) = f £„(r(r)) exp(-r)dr + (1 - re,. - tr c )B„(T e ) exp(-ri)+ 

re c exp(-n) /“ B„(T(r)) Htp(-r)dr+tr t (l-re<) «p(-r,) [” B„ (T(r)) exp(-r)dr 

(3.73) 

3.3.3 Vertical Distributions of Opacity Sources 

In the previous two sections, we constrained two of the three parameters of the 
RTM to a large degree. This leaves the vertical distribution of opacity sources 
as the main variable in the model. The gaseous constituents are well mixed in 
the deep atmospheres of the giant planets. Some of the gaseous constituents will 
be destroyed by photolysis or condensation at higher altitudes. We estimate the 
vertical distribution of the cloud forming constituents NH S , H 2 S, and H 2 0 using the 
thermochemical model described in Section 3.1. The actual distributions are more 
complex; they exhibit latitudinal variations due to upwelling and downwelling, and 
some of the gaseous constituents might be super-saturated (see, e.g., de Pater et 
al., 1989). We will use a single distribution to model the disk-averaged brightness 
temperatures of the giant planets. Our goal is to find distributions of opacity 
sources which provide good fits to the observed spectra of giant planets. 


3.4 Modeling Results 

In the previous sections, we developed an elaborate model for computing the 
brightness temperatures of the giant planets. We now adjust the parameters of the 
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model in order to obtain a good fit between the modeled and observed brightness 
temperatures. A list of reliable observations has been compiled for each of the 
planets. We examine the millimeter observations and their calibration procedures 
in detail. The most detailed analysis is carried out for Jupiter. Of the four gi- 
ant planets, Jupiter has perhaps the most interesting millimeter-wave spectrum 
because the emission originates from pressures near the putative NH 4 SH cloud. 

The temperature-pressure profile, T(P) for each of the Jovian planets is shown 
in Figure 3.8. The temperature-pressure profiles of the four planets are similar in 
that each planet has a well defined troposphere in which the temperature increases 
with increasing pressure. The temperatures of the planets generally decrease as 
the distance from the sun increases. The temperature-pressure profiles for Uranus 
and Neptune are almost identical even though Neptune is much further from the 
sun than Uranus. Neptune, like Jupiter and Saturn, possesses a large internal 
heat source. Uranus has a much smaller internal heat source. Some of the model 
parameters are listed in Table 3.2. The deepest temperature and pressure probed 
by the Voyager spacecraft are listed in the table. 

The modeling work presented here extends the work of others who have de- 
veloped similar models to carry out detailed analyses of the centimeter wavelength 
emission from the giant planets (see references in following sections). The goal here 
is to focus on the millimeter region of the radio spectrum and to draw comparisons 
between the millimeter- wave spectra of the four giant planets. 

3.4.1 The radio spectrum of Jupiter 

Jupiter’s radio emission is composed of both thermal (quasi-thermal) and non- 
thermal components. The thermal component is emission from the planet which 
exhibits the characteristics of a blackbody radiating at its physical temperature. 
The non-thermal component is cyclotron and synchrotron radiation from energetic 
electrons moving within the magnetosphere of a planet. The non-thermal com- 
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Temperature- Pressure Profiles of Planets 



Figure 3.8: Temperature-pressure profiles of the Jovian planets; Solid line: Jupiter, 
dashed line: Saturn, dot-dashed: Uranus, and dotted line: Neptune (indistinguish- 
able from Uranus). 
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Table 3.2: Model parameters 



Jupiter 

Saturn 

Uranus 

Neptune 

Xh 3 

0.90 

0.963 

0.833 

0.833 

g (cm/s J ) 

2417 

982 

850 

1130 

T (K) 

165° 

107.5* 

100.9* 

78* 

P (bar) 

1 

0.5 

2.3 

1.2 


‘Lindal et a/.(1981), *Lindal et o/.(l985),*Lindal et o/.(1987), ‘'Tyler et al.{ 1989) 


ponent is dominant at wavelengths longer than 7 cm and negligible at millime- 
ter wavelengths. The non-thermal component can be separated from the thermal 
component so that the observed thermal emission may be compared with synthetic 
thermal emission spectra. 

We have compiled a list of Jupiter’s published millimeter observations to 
be used as a basis for comparison with synthetic spectra. Table 3.3 lists the 
reliable observations at frequencies between 36 and 300 GHz (1-8.35 mm) along 
with corresponding calibration sources and spectral bandwidths. Klein and Gulkis 
(1978) have normalized all of the observations between 14.5 and 36 GHz (8.35-20.7 
mm) to a common flux scale. At longer wavelengths, Berge and Gulkis (1976) have 
given a survey in which the observations have been normalized to a common flux 
scale whenever possible. They have also removed the non-thermal component of 
Jupiter’s flux. 

Vertical distributions of NH 3 , HjS, H 2 O, and PH 3 are shown in Figure 3.9. 
The distributions of NH 3 , HjS, and H 2 0 are estimated with the thermochemical 



Table 3.3: List of reliable millimeter observations of Jupiter 


A 

(mm) 

V 

(GHz) 

Tb 

(K) 

Cal 

Ai/ 

(GHz) 

Reference 

1.0 

300 

168±8 

Melts 

102,229 

Werner et al. (1978) 

1.08 

279 

169.9±5.1 

Mars 

75 

Griffin et a/.(1986) 

1.32 

227 

170.9±3.9 

Mars 

70 

Griffin et al. (1986) 

1.32 

227 

165±8 

planets 

39 

Ulich et a/. (1984) 

1.4 

214 

148±16° 

planets 

275 

Rather et al. (1975) 

1.40 

214 

168±‘ll 

Mars 

210 

Courtin et of. (1977) 

2.00 

150 

173.3±1.1 

Mars 

50 

Griffin et al. (1986) 

2.13 

141 

167±12“ 

abs 

1-2 

Ulich (1974, 1981) 

2.14 

140 

178±13“ 

abs 

1-2 

Cogdell et al. (1975) 

3.09 

97 

174±10° 

abs 

0.1-0.2 

Ulich et al.{ 1973) 

3.33 

90 

172.5±1.4 

Mars 

1-2 

Griffin et al.( 1986) 

3.48 

86 

179.4±4.7 

abs 

0.03-0.5 

Ulich et of. (1980) 

3.53 

85 

166±6 B 

DR21 

1-2 

Ulich (1974) 


“Brightness temperatures as given in Berge and Gulkis (1976) 
‘Recalculated with beam correction factor in Ulich et o/.(1980) 
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model described in Section 3.1. The NH 3 vertical distributions assume a deep 
mixing ratio of 2.5 X 10 -4 . This value corresponds to the upper limit derived by 
Bjoraker tt al. (1986a) at 5 fim. Distributions of H 2 S assume subcloud mixing 
ratios equal to the solar abundance (distribution 1: XH a s = 2.9 x 10 -5 , Anders and 
Grevesse, 1989) and approximately seven times the solar abundance (distribution 
2: XH a s = 2.2 x 10 -4 ). Since the H 2 S abundance is less than the NH 3 abundance, 
H 2 S is rapidly depleted near 2 bar due to the reaction with NH S to form NH4SH. 
The upper limit for the deep mixing ratio of H 2 S is the deep mixing ratio of NH 3 . 
If the mixing ratio of H 2 S were greater than the m ixin g ratio of ammonia, the 
ammonia would be rapidly depleted at 2 bar. This would be inconsistent with 
microwave observations (Lindal tt al ., 1981 and de Pater and Massie, 1985). The 
H 2 S distributions are compatible with the upper limit derived by Larson tt al. 
(1984) at 2.7/xm. For subcloud H 2 S mixing ratios less than or equal to the solar 
abundance, NH 3 is not significantly depleted (as in NH 3 distribution 1). However, 
if the H 2 S mixing ratio is increased to 2.2 x 10 -5 (as in H 2 S distribution 2), NH 3 
becomes depleted near 2 bar (as in NH 3 distribution 2). The formation of an 
NH 3 ice cloud further depletes the ammonia near 0.75 bar. The H 2 0 distribution 
assumes a subcloud mixing ratio equal to the solar abundance (1.5 x 10~ s ) and is 
depleted by the formation of an H 2 0 cloud. The distribution of PH 3 is the same as 
that used by Lellouch tt al. (1984c) which is based on the inferred profile of Kunde 
tt al. (1982). This distribution assumes a deep PH 3 mixing ratio xph* = 6 x 10~ 7 . 
Phosphine is depleted by photolysis above 1 bar. 

The distributions and cloud bulk densities shown here may not represent the 
actual distributions in Jupiter’s atmosphere, but are used as a basis to test the 
effects of the various absorbers on the computed Jovian spectrum. For example, 
the actual H 2 0 distribution in Jupiter’s atmosphere may be depleted by a factor 
of 100 (Bjoraker tt al. t 1986b). The H 2 0 distributions and bulk densities shown 
above are used to illustrate the maximum possible effect of H 2 0 (gas and clouds) 



VOLUME MIXING RATIO/CLOUD BULK DENSITY g/ml 
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Figure 3.9: Vertical distributions of NH S , H,S, PH S , and H,0 in Jupiter’s at- 
mosphere. XNE t = 2.5 x lO" 4 , XHjO = 15 x 10 -* (solar), Xpe. = 6 x 10 " 7 ; 
Distribution Is Xb,s = 2.9 x 10~ l (solar), Distribution 2: Xh,s = 2.2 x 10 " 4 . 
Cloud bulk densities from thermochemical model using H,S distribution 2. 
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on Jupiter’s emission. We do not consider the latitudinal variations in the cm 
wavelength emission which have been observed with the VLA by de Pater (1986). 
Instead, we use a single vertical distribution to represent an average over the disk 
of the planet. 

Model sensitivity to NHs and H2O 

In Figure 3.10, the Jupiter’s observed spectrum is shown along with computed 
emission spectra. The synthetic emission spectra include only NH S , H 2 O, and 
pressure-induced opacity and use the two NHj distributions shown in Figure 3.9. 
The error bars for observations using Mars as the primary calibrator include a 
10% uncertainty for the calibration in addition to the systematic uncertainties 
reported in Table 3.2. Synthetic spectra using both NH S distributions provide a 
good fit to observations near 1.3 cm and at longer wavelengths. However, only 
NHj distribution 1 provides a good fit to the millimeter-wave observations. The 
influence of the 183 GHz HjO line on the computed emission is small for NHs 
distribution 1. 

Model sensitivity to H2S and PH3 

The effect of adding H 2 S opacity to Jupiter’s computed spectrum for the two NHs 
distributions is shown in Figure 3.11 (dashed line: NH S , H 2 0, and pressure-induced 
opacity only, solid line: H 2 S opacity added). We have also modeled the effect of 
the phosphine (PH S ) rotational line at 267 GHz. As seen in Figure 3.12, the effect 
of PH S opacity is limited to frequencies near its transition at 267 GHz. 

Model sensitivity to clouds 

Previous attempts to model the effects of cloud opacity have only included the po- 
tential contribution from the NHs ice cloud. For example, Griffin et al. (1986) have 
modeled the potential effect of the NHs cloud on Jupiter’s millimeter and submil- 
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Figure 3.10: Observed spectrum of Jupiter. Error bars include systematic uncer- 
tainty (solid line) and 10% uncertainty for observations using Mars as calibrator 
(dashed line). Computed emission using NH», H a O, and pressure-induced opacity 
only with vertical distributions in Figure 3.8. 




BRIGHTNESS TEMPERATURE Kelvins 
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Figure 3.11: Observed spectrum of Jupiter with computed emission using NH S and 
HjS distributions 1 and 2 from Figure 3.9. Dashed line: NHs and HjO opacity 
only; Solid line: NH S , H a O, PH S , and HjS opacity. 
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limeter emission using a multiple scattering solution for various effective radii and 
particle- to-gas scale heights. Their model assumed that the imaginary part of the 
index of refraction (k NHt ) decreased exponentially with frequency from the mea- 
sured value at 50 cm" 1 (A = 200 fim). The difference in brightness temperature 
resulting from cloud scattering model at 1 mm was less than 5%. Similarly, Bezard 
et a/. (1983) have modeled the multiple scattering effects of the NH S cloud by as- 
suming that kffn» was the same at 1 mm as the measured value at 200 /xm. They 
also found that the contribution from the NH S cloud at millimeter wavelengths 

was negligible. 

Figure 3.12 shows normalized weighting functions at 1 mm, 1 cm, and 10 cm 
using model 2 above with and without cloud opacity. The 1 cm weighting function 
is very sensitive to gaseous NH S since it is near the 1.3 cm inversion wavelength. 
The 1 cm weighting function peaks near 0.5 bar, close to the NH S saturated vapor 
pressure. At longer wavelengths, NH S is more transparent, and the weighting 
function peaks much deeper in the atmosphere. For example, the 10 cm weighting 
function peaks near 6 bar. At millimeter wavelengths, the weighting functions 
begin to increase near 0.8 bar due to NH 3 opacity. The 1.4 mm weighting function 
has a second peak near 2 bar. This second peak is due to H 2 S absorption since 

1.4 mm is near a strong H 2 S absorption line. 

The weighting functions show that most of the emission at short millime- 
ter wavelengths originates from pressures between 1 and 3 bars. Therefore, the 
potential contributions from the putative NH 4 SH cloud should also be included. 
We have computed a weighting function which includes cloud opacity using Equa- 
tion 3.56. This weighting function does not include the effects of cloud scattering. 
The 1 mm weighting function with cloud opacity has several interesting features. 
The main peak near 0.8 bar is due to gaseous NH s (also present in the cloud free 
weighting function). The first small peak near 0.5 bar is due to opacity from the 
NH S ice cloud. The second much larger peak is due to opacity from the NH 4 SH 


Weighting Functions (arb units) 
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Weighting functions for Jupiter 



PRESSURE bar 


Figure 3.12: Weighting functions for Jupiter at 1 mm (with and without cloud 
opacity, 1.4 mm (with HjS opacity), 1 cm, and 10 cm. 
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cloud near 2 bar. 

We have computed the cloud reflectivity (re tf ) and transmissivity (tr e ) of the 
NH S and NH 4 SH clouds using the approach outlined in Section 3.3.2 for scat- 
tering in the forward and backward directions. We compute the reflectivity and 
transmissivity for several different values of k (the imaginary or absorbing part 
of the complex index of refraction, n = n - jfc), fractional bulk density M/M oy 
where M 0 are the bulk densities shown in Figure 3.9, and mode radius (r e ). The 
bulk densities shown in Figure 3.9 were computed by assuming that the subcloud 
H 2 S abundance was enhanced relative to the solar value. If the H 2 S abundance 
is closer to the solar abundance, the NH 4 SH cloud density will be lower and the 
NHs cloud density will be higher. Therefore, we have computed the reflectivities 
and transmissivities over a range of bulk densities. The computed reflectivities, 
transmissivities, and corresponding zenith optical depths, t = -ln(t e ), are listed 
in Table 3.4. The transmissivity includes absorption from the clouds as well as 
the gases within the cloud. We have also determined the corresponding reduction 
in Jupiter’s brightness temperature resulting from the effects of reflection and ab- 
sorption by clouds using the two NH 3 distributions shown in Figure 3.9 [AT, and 
A? 2 , respectively) . The clouds have a larger effect when less opacity is present 

above the clouds as in NHj distribution 2. 

Both the opacity and reflectance of the NH, ice cloud are small. We also 
find that its influence on Jupiter’s spectrum is small. However, the opacity and 
reflectance of the NH 4 SH ice cloud is significant in some cases. For r e > 40 /im, 
the cloud is highly absorbing regardless of the value of k NHiS n- For r e = 10 M m, 
the cloud is still highly absorbing for k N H t sH > 001 • The NH < SH cloud is also 
opaque when the cloud bulk density is decreased by a factor of 20. 

Figure 3.13 shows synthetic emission spectra with and without a highly ab- 
sorbing cloud at 2 bar using NH 3 distributions 1 and 2. The absorption is com- 
puting according to Equation 3.56 (scattering is not included) with k = 0.05 and 
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Table 3.4: Potential cloud reflectivities (re e ), transmissivities (tr e ), zenith optical 
depth (r), and reduction in Jupiter’s brightness temperature using NHj distribu- 
tions 1 and 2 from Figure 3.9 (ATi and AT]) at 1 mm for different values of k, 
n = n — jk, mode radius r e , and fractional bulk density M/M 0 where M 0 are the 
bulk densities shown in Figure 3.9 


Cloud k 

M/M 0 

r 0 (nm) 

tr« 

T 

re e 

AT a 

AT] 

NHj 0.001 

1 

40 

0.99 

0.01 

0.0008 

1 

3 

(n = 1.3) 0.01 

1 

40 

0.98 

0.02 

0.005 

3 

5 

0.05 

1 

40 

0.90 

0.11 

0.008 

8 

12 

0.001 

1 

100 

0.99 

0.01 

0.002 

1 

3 

0.01 

1 

100 

0.97 

0.03 

0.01 

5 

7 

0.05 

1 

100 

0.89 

0.12 

0.02 

12 

16 

0.001 

2 

40 

0.99 

0.01 

0.001 

1 

3 

0.01 

2 

40 

0.95 

0.05 

0.01 

6 

9 

0.05 

2 

40 

0.81 

0.21 

0.02 

16 

22 

0.001 

2 

100 

0.99 

0.01 

0.01 

4 

6 

0.01 

2 

100 

0.93 

0.07 

0.03 

13 

17 

0.05 

2 

100 

0.79 

0.24 

0.04 

22 

29 

NH 4 SH 0.001 

1 

10 

0.16 

1.8 

0.06 

5 

17 

(n = 1.74) 0.01 

1 

10 

0.01 

4.6 

0.03 

4 

14 

0.05 

1 

10 

5xl0 -8 

17 

0.01 

4 

11 

0.001 

1 

40 

2xl0 -4 

8.5 

0.36 

4 

55 

0.01 

1 

40 

4 x 10 -8 

17 

0.28 

11 

45 

0.05 

1 

40 

4 x 10 -16 

35 

0.19 

10 

33 

0.001 

1 

100 

lxKT 6 

14 

0.60 

8 

85 

0.01 

1 

100 

1 x 10 -13 

30 

0.46 

17 

68 

0.05 

1 

100 

9xlO“ M 

62 

0.26 

14 

43 

0.001 

0.05 

10 

0.30 

1.2 

0.005 

2 

8 

0.01 

0.05 

10 

0.26 

1.3 

0.004 

2 

8 

0.05 

0.05 

10 

0.14 

2.0 

0.003 

3 

9 

0.001 

0.05 

40 

0.15 

1.9 

0.12 

6 

26 

0.01 

0.05 

40 

0.13 

2.0 

0.11 

5 

25 

0.05 

0.05 

40 

0.05 

2.7 

0.08 

5 

20 

0.001 

0.05 

100 

0.07 

2.7 

0.21 

9 

34 

0.01 

0.05 

100 

0.05 

3.0 

0.19 

8 

35 

0.05 

0.05 

100 

0.01 

4.6 

0.13 

6 

26 



81 


the bulk densities equal to those shown in Figure 3.9 in order to simulate a highly 
absorbing NH 4 SH cloud. At 1 mm, the brightness temperature using both distri- 
butions with a highly absorbing cloud provide good fits to the observations. The 
clouds also decrease the effect of HjS shown with a solid line, since most of the 
HjS gas should exist just beneath the NH^SH cloud. Since neither cloud nor NHj 
absorption has distinctive spectral features, cloud opacity cannot be distinguished 
from the effects of NH S opacity. However, NHj distribution 1 still provides a bet- 
ter fit to the spectrum at the longer millimeter wavelengths than distribution 2. 
This suggests that NH S is most likely the dominant source of opacity on Jupiter 
at millimeter wavelengths. This also suggests that NH, is probably not signifi- 
cantly depleted by the reaction with H 2 S to form solid NH 4 SH and consequently 

XH a S < Xnh 5 - 

3.4.2 The radio spectrum of Saturn 

A list of Saturn’s radio observations has been compiled by Klein ct al. (1978). The 
observations in this survey have been normalized to a common flux scale. Correc- 
tions for the size of the oblate planet were also applied to all of the observations. 
The influence of the rings was removed so that the brightness temperatures corre- 
spond to the temperatures which would have been observed in the absence of the 
rings. Additional observations have been compiled by Grossman (1990). 

Detailed analyses of the centimeter spectrum of Saturn have been carried out 
by Grossman (1990) and Briggs and Sackett (1989). These authors have investi- 
gated several possible distributions of NH S , HjS, and H 2 0 in Saturn’s atmosphere. 
Their best fits to centimeter wavelengths observations used enriched NH S and H 2 S 
mixing ratios (relative to the sun). As on Jupiter, the overall mixing ratio of H 2 S 
is less than that of NH S , so that H 2 S is depleted in the upper atmosphere. The 
observed latitudinal variation at centimeter wavelengths is characterized by a hot 
band at northern latitudes (Grossman, 1990). This is suggestive of a large scale 



BRIGHTNESS TEMPERATURE Kelvins 


82 


Jupiter’s spectrum with cloud opacity 



WAVELENGTH cm 


Figure 3.13: Observed spectrum of Jupiter with computed emission using NH S 
distributions 1 and 2 with and without cloud opacity included. Dashed line: NH S 
and HjO opacity only; dot-dash line: highly absorbing cloud; solid line: HjS 
opacity added. 



83 


circulation pattern. 

Figure 3.14 shows three vertical distributions of NH S and cloud bulk densities. 
Figure 3.15 shows Saturn’s observed spectrum with three synthetic spectra using 
the vertical distributions in Figure 3.14. Figure 3.16 shows normalized weighting 
functions using NH S distribution 3 at 1 mm, 1 cm, and 10 cm. Again, the weighting 
function at 1 cm is sensitive to the presence of NH S gas and peaks near the NH 3 
saturation pressure. The weighting function at 1 mm peaks slightly deeper in 
the atmosphere. The weighting function at 10 cm peaks still deeper. The sharp 
peak in the 10 cm weighting function near 7 bar is due to the increase in the NH S 
abundance at that pressure (see Figure 3.14) where the NH4SH cloud forms. 

Since the millimeter-wave emission originates from near the same altitudes 
as the NH S cloud, Saturn’s cloud opacity should be very small. Moreover, the 
millimeter-wave emission originates from altitudes where H 2 S and H 2 0 are likely 
to be severely depleted. Therefore, the only opacity sources are NH S and pressure- 
induced absorption. Since the opacity from NH S is significantly larger than that 
from pressure-induced absorption, the millimeter-wave spectrum of Saturn is pri- 
marily sensitive to the vertical distribution of NHj. 


3.4.3 The radio spectrum of Uranus 

Gulkis and de Pater (1984) have reviewed the existing millimeter and centimeter 
observations of Uranus and have normalized the observations to a common radius 
and flux scale whenever possible. Detailed analyses of the microwave spectrum of 
Uranus have been carried out by Hofstadter and Muhleman (1988), Hofstadter et 
at. (1990) and de Pater et al. (1989, 1991). Both temporal and spatial variations 
have been observed in the Uranus atmosphere (see, e.g., , Hofstadter and Muhle- 
man, 1988 and Gulkis and de Pater, 1984). The spatial variation is characterized 
by a large pole-equator gradient in which the brightness temperature is warmer 
at the pole than at the equator. The observed brightness temperatures of Uranus 
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PRESSURE bar 


Figure 3.14: Vertical NH S distributions and cloud bulk densities in Saturn’s at- 
mosphere; 1-Dashed line: xni r» = 4.5 x 10 -4 , Xh 3 s = 1 x 10 -4 ; 2-Dot-dashed line: 
Xnh, = 1.5 x 10 -4 , XFjS = 2.9 x 10 -1 (solar); 3- Solid line: Xnh, = 4 x 10 -4 , 
Xe 3 s = 3.6 x 10 -4 ; For cloud bulk densities, Xnh» = 4 x 10 -4 , xs a s = 3.6 x 10 -4 , 
and xb»o = 1.5 x 10 -s . 





PRESSURE (bar) 


Figure 3.16: Weighting functions at 1 mm, 1 cm, and 10 cm for Saturn using NH< 
distribution 3. 
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increased between 1965 and 1973. This was the time when the its pole, inclined 
at an angle of almost 90 degrees, came into view. The temporal variation can be 
explained to first order by a change in the viewing aspect of the planet (see t.g., 
de Pater tt al., 1989). 

The observed radio spectrum of Uranus is shown in Figure 3.17 along with 
several synthetic spectra. The spectrum at short millimeter wavelengths is domi- 
nated by pressure-induced absorption. The vertical distributions and mixing ratios 
of NH S , H 2 S, and HjO have virtually no effect on Uranus’ millimeter spectrum. 
Therefore, our computed spectra should in effect be a good test of the assumed 
temperature-pressure profile. At centimeter wavelengths, we show only a few of 
the computed spectra which have been previously explored (see, t.g., de Pater tt 
al., 1989, 1991). We show the effect of pressure-broadened HjS on Uranus’ cen- 
timeter wavelength emission. The synthetic spectra with xh 7 s > Xnh» provide a 
better fit to the observed spectrum than those with xh 9 s < Xnh 9 - 

Weighting functions at 1 mm, 3 mm, 1 cm, and 10 cm are shown in Fig- 
ure 3.18. The 1 cm weighting function has a sharp peak due to NH S opacity. This 
peak occurs at higher pressures than on Jupiter and Saturn due to the colder 
temperatures. The weighting function at 10 cm is similar to that of Saturn. The 
second peak near 20 bar is due to an increase in NH S abundance near the saturation 
pressure of NH 4 SH. In contrast to Jupiter and Saturn, the 1 mm weighting peaks 
at higher altitudes (lower pressures) than the 1 cm weighting function. This is 
evidence of the large effect of pressure-induced opacity at millimeter wavelengths. 

3.4.4 The radio spectrum of Neptune 

de Pater and Richmond (1990) have compiled a list of Neptune’s radio observa- 
tions. The microwave observations in this list were normalized to the common 
flux scale of Baars tt al. (1977). The millimeter observations which were cali- 
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Uranus' Radio Spectrum 


no H„S 


WAVELENGTH cm 


Figure 3.17: Observed spectrum of Uranus with computed emission; la, Solid line: 
Xnh* = 1*5 x 10 -s , Xh»s = 1-Ox 10~ 4 (no H-.S opacity);lb, Dashed line: H 3 S opacity 
added; 2, Dash-dot line: xnh, = 4*0 x 10 _B , xh*s = 3.3 x 10 - *; xh,o = 1.5 x 10 -8 
for all spectra. 
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brated with Mars or Jupiter were renormalized following the method described in 
Griffin et al. (1986). 

Neptune’s millimeter spectrum is quite similar to that of Uranus. At short 
millimeter wavelengths, the opacity is primarily due to pressure-induced absorp- 
tion. Again, the millimeter spectrum is relatively insensitive to the mixing ratios of 
the volatile constituents. Romani et al. (1989) and de Pater et al. (1991) have an- 
alyzed Neptune’s centimeter spectrum in great detail. At centimeter wavelengths, 
the spectrum is dominated by NH S , liquid H 2 0, and possibly HjS opacity. We 
show Neptune’s observed spectrum and a sample of computed spectra which are 
similar to those previously explored by Romani et al. (1989) and de Pater et al. 
(1991) in Figure 3.19. 


3.5 Conclusions and Comparisons 

The weighting functions at 1 mm peak near 1 bar for all four planets. Table 3.5 
lists the magnitude of the absorption coefficients for the four Jovian planets at 1 
bar and at 1 mm. At millimeter wavelengths, Jupiter and Saturn exhibit char- 
acteristics of NH 3 opacity. The millimeter-wave spectra of Uranus and Neptune 
are dominated by pressure-induced absorption and not NH 3 opacity. They do 
not provide any information about NH 3 or HjS abundances or distributions. All 
four planets exhibit higher brightness temperatures at the longer centimeter wave- 
lengths which is primarily due to the inversion spectrum of NH 3 . 

Using our new formalism for the Ben-Reuven line shape in the radiative 
transfer model, we found that Jupiter’s observed millimeter-wave spectra can ad- 
equately be explained with NH 3 as the primary opacity source. Good fits to the 
observed spectrum are achieved using vertical distributions for NH 3 derived from 
equilibrium condensation models with Xsn t = 2.5 x 10 -4 (solar abundance en- 
hanced by a factor of 1.3) and Xb 3 s < 2.9 x 10~ 5 (solar abundance). 
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Figure 3.10: Observed spectrum of Neptune with computed emission; la, solid 
line: xnh. = 15 X 1CT 4 , X h,s = 3.3 x 10* 4 (no H,S opacity); lb, Dashed line: 
H,S opacity added; 2, Dash-dot line: xnb, = 4 0 x 10 - ', Xir.s = 3.3 x 10" B ; 
XH,o = 1-5 x 10 -8 for all spectra. 
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Table 3.5: Magnitude of absorption coefficients (in cm *) at 1 bar, 1 mm 



Jupiter 

Saturn 

Uranus 

Neptune 

&H 2 

lx nr 8 

3xl0 -8 

2xl0 -7 

2xl0 -7 

a NH» 

3xl0 -7 

lxlO -7 

lxlO -13 

lxlO" 13 

<*h 2 s 

8xl0 -11 

lx icr 10 

lxlO -17 

lxlO -27 

<*B,0 

2xl0 -u 

4xl0 -11 

2x10" 26 

2xl0~ 16 

<*PH» 

4xlO -10 

- 

- 
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The influence of gaseous HjO on Jupiter’s spectrum is small. Detection 
of the 183 GHz HjO line would be extremely difficult with ground-based radio 
telescopes due to blanking by the terrestrial HjO line. Cloud condensates may be 
providing opacity. However, due to large uncertainties in the dielectric properties 
of the condensates and cloud bulk densities, no firm conclusions regarding cloud 
opacities can be drawn at this time. Gaseous hydrogen sulfide may be providing 
additional opacity on Jupiter. Observations with greater spectral resolution are 
needed in order to detect the potential effect of HjS opacity on Jupiter’s millimeter 
spectrum and thus distinguish between the effects of HjS and NH 3 opacity. We 
discuss the approach and results of such an observation in the following chapter. 
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CHAPTER 4 

Dual Wavelength Observation of Jupiter 

at 1.4 mm 


Sulfur is the tenth most abundant element in the solar system (Anders and Grevesse, 
1989). Although sulfur has been detected on Jupiter’s volcanic moon Io, it remains 
mysteriously absent from the list of elements detected in the atmospheres of the 
Jovian planets. Ground-based observation with radio telescopes provides one of 
the few means capable of probing beneath the dense cloud layers of the Jovian 
planets where sulfur in the form of gaseous hydrogen sulfide (H 2 S) may be found. 
Jupiter presents the best target for the H 2 S search due to its cloud structure and 
extreme brightness at millimeter wavelengths. Determination of the sulfur abun- 
dance in the giant planet atmospheres would provide valuable information needed 
to develop and refine models of the origin and evolution of the planets and so- 
lar system - The search for sulfur is also of interest relative to the long-standing 
problem of the unknown composition of colored cloud materials in Jupiter. Be- 
cause amm onium polysulfides and other sulfur compounds are leading candidates 
for this material (West et a/., 1986), detection of any sulfur containing molecule 
would be of great value. 

Both theory and observation suggest that nearly all of the sulfur is locked up 
in the clouds of Jupiter’s troposphere. Thermochemical models predict that H 2 S 
will combine with NH S to form an ammonium hydrosulfide (NH 4 SH) ice cloud near 
the 2 bar pressure level on Jupiter (see Section 3.4.1). If the abundance of H 2 S is 



94 


less than the abundance of NH 3 at subcloud altitudes in Jupiter’s atmosphere, the 
reaction to produce NH4SH should rapidly deplete H2S. Any HjS which escaped 
destruction in the reaction to form NH4SH should be destroyed by photolysis near 
1 bar (Bezard tt al., 1983). Infrared observations of Jupiter at 2.7 nm by Larson 
tt al. (1984) provide a strict upper limit on the abundance of H2S at pressures 
between 0.7 and 1.2 bar. This limit reinforces the hypothesis that nearly all of the 
H2S is depleted by the formation of an NH4SH cloud. This leaves millimeter-wave 
observations as one of the only means of probing deep enough in the atmosphere 
to detect H2S. 

Bezard tt al. (1983) suggested that the rotational lines of gaseous hydrogen 
sulfide (H2S) in Jupiter’s troposphere might provide excess opacity in Jupiter’s 
millimeter spectrum. Millimeter-wave observations probe altitudes in the Jovian 
atmosphere at and below the altitude of the putative NH4SH cloud base, where 
gaseous H 2 S may exist. As shown in Figure 3.11, three strong rotational lines of 
H2S could potentially affect the emission from Jupiter at millimeter wavelengths. 
The center frequencies of these lines are 168.8, 216.7 and 300.5 GHz (1.78, 1.38, and 
1.00 mm). The hydrogen sulfide (H 2 S) transition at 216.7 GHz provides the best 
opportunity to search for H2S because it appears in relatively transparent regions 
of the atmospheres of Jupiter and the earth. In Jupiter’s atmosphere, opacity from 
NH 3 obscures the effects of potential H 2 S lines at wavelengths shorter than 1 mm. 
Interference from the 183 GHz water vapor line in the earth’s atmosphere would 
complicate an observation of the 168.8 GHz H2S transition. 

Detecting the 1.4 mm (216.7 GHz) H2S line with existing instrumentation 
is difficult. This line is pressure-broadened by hydrogen (H s ) and helium (He), 
spreading its total width over 30-40 GHz. Most of the existing observations of 
Jupiter at millimeter wavelengths have been made with broadband filters which 
have pass bands of approximately 70 GHz. The filters effectively average the 
emission detected in the pass band, smearing any smaller scale spectral features 


95 


such as pressure-broadened H,S lines. Other receivers used in millimeter and 
submillimeter astronomy are designed in order to detect narrow emission lines 
( A „/„ ~ 10 -*) in the interstellar medium. These receivers typically have pass 
bands of less than 1 GHz which would be too narrow to detect pressure-broadened 

H 2 S lines. 

We have developed an approach in which a high resolution receiver can be 
used in order to search for broad spectral features, such as the pressure-broadened 
H 2 S lines on Jupiter. This technique involves observing the differential emission 
between two or more wavelengths in order to detect H 2 S or place limits on its 
abundance. We have used this approach to observe Jupiter with the 10.4 m Cal- 
tech Submillimeter Observatory (CSO) 1 at Mauna Kea, Hawaii. In addition, we 
have used Mars as a calibration source in order to derive a reliable brightness 

temperature for Jupiter at 1.4 mm. 

4.1 Sensitivity Calculation 

Before any observations were undertaken, we had to confirm that the receiver 
would be sensitive to the change in Jupiter’s emission due to H 2 S absorption. It 
is useful to first define several quantities used in this type of calculation. The 
Rayleigh- Jeans approximation is commonly used at radio frequencies m order to 
simplify Planck’s law (Equation 3.10): 

B = 2i SkT/c 1 = 2kT/X 7 . ( 4>1 ) 

The Rayleigh-Jeans approximation is valid for hu < kT. The approximation pro- 
duces errors of approximately 10% at millimeter wavelengths. Using the Rayleigh- 
Jeans approximation, the total spectral power density, 5, radiated from a spherical 
' iThe CSO is operated by the California Institute of Technology under funding from the National 
Science Foundation, Contract #AST-88-15132. 
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blackbody of radius r at a distance d from the source is 

S = 27rfcr(r/d)VA 3 , (4.2) 

(Gulkis, 1987), where S is called the flux density. The flux density is in units of 
power per unit surface area per unit frequency. A commonly used unit for flux 
density is the Janksy (Jy) which is defined as 10” W m _J Hz -1 . The change in 
antenna temperature, AT, (in Kelvins) due to a certain flux density, 5, is 

AT = (4.3) 

where A,® is the effective area of the antenna. The effective area is 

R 3 , (4.4) 

where a is the efficiency of the dish, and R is the radius of the antenna. 

We first compute the change in antenna temperature (A7# a s) resulting from 
the expected 2 K difference in Jupiter’s brightness temperature due to H 2 S ab- 
sorption using Equations 4.2 -4.4. For Jupiter, r ~ 71.6 x 10 s km and d c* 4.2 AU 
(astronomical units) or 6.28 x 10 8 km. The diameter of the CSO is 10.4 m and an 
efficiency of 30% is assumed. We computed the difference in emission assuming a 
temperature of 180 K at 216 GHz and a temperature of 182 K at 200 GHz. After 
taking into account the differential flux due to the difference in wavelength, ATjj 3 s 
is 0.76 K. 

We now examine the sensitivity of the receiver in order to determine whether 
or not it is capable of measuring ATg t s- The sensitivity equation of an ideal 
balanced Dicke radiometer is 

AT ™ = 2 ( 4 5 ) 

where ATnn* is called the rms noise power, T$ is the system temperature (a measure 
of the noise power from the receiver), t is the integration time in seconds, and du 
is the band width of the receiver in Hz. The full receiver band width of the CSO is 
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500 MHz. The receiver operates between 200 and 260 GHz. Although the CSO is 
equipped to make observations with much greater resolution, the full band width 
of the receiver will be utilized in order to achieve the necessary sensitivity. 

The noise temperature of the CSO (for double side band) is reported to be 
100 K at 200 GHz and the system sensitivity approximately 500 mJy for a one 
second integration time. Evaluating Equation 4.5 gives a AT^, < 0.1 K. Since 
ATnn, < AThjS, we are not limited by the sensitivity of the CSO receiver. 


4.2 Instrumentation and Procedure 

Figure 4.2 shows a block diagram of the CSO facility. The CSO has two available 
receivers. The 345 GHz receiver is located at the Cassegrain mount. The 230 
GHz receiver is located at the sidecab mount in a separate room. A mirror is 
used to optically transmit the signal from the antenna to the 230 GHz (sidecab) 
receiver. Our observations were made with the 230 GHz receiver. The receiver is 
a superconducting insulator superconducting (SIS) quasiparticle tunnel junction 
mixer. It is contained in a liquid helium dewar at a temperature of about 4 K. A low 
noise amplifier is maintained at a temperature of about 12 K. The receiver operates 
in a double side band mode with a band separation of 2.8 GHz (IF band 1150- 
1650 MHz). A detailed drawing of the 230 GHz receiver is shown in Figure 4.3. 
The receiver’s LO (local oscillator) is generated by doubling a signal from a Gunn 
oscillator. The oscillator is phase locked to a microwave frequency counter. The LO 
signal is optically coupled to the receiver through a horn antennas and focusing 
lenses. The CSO is equipped with two available Acousto-Optic Spectrometers 
(AOS). Each AOS has 1024 channels with total bandwidths of 500 or 50 MHz. The 
backend computer automatically controls attenuator settings in order to prevent 
saturation of the AOS, 

The receiver requires careful tuning at each frequency. Tuning the receiver 
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Figure 4.1: Block diagram of the Caltech Submillimeter Observatory (CSO). 
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Signal from Dish (Plonet) 


Figure 4.2: Block diagram of the CSO 230 GHz receiver 
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involves the manual adjustment of several back shorts and attenuators in order to 
optimize the receiver temperature (see Figure 4.3). Because the two frequencies we 
chose were separated by approximately 15 GHz, it was not possible to use a single 
local (Gunn) oscillator (LO). Changing oscillators increased the time required to 
tune the receiver. We required between 45 minutes and 2 hours to achieve adequate 
receiver temperatures. 

The receiver temperature can be computed by measuring the voltages at the 
receiver with a hot load (piece of absorber at ambient temperature T k ^ 273 K) 
and cold load (piece of absorber dipped in liquid nitrogen at a temperature T e ~ 
77 K). The hot and cold load voltages (V* and V e , respectively) are related to the 
receiver temperature (Trx) by 


V* = G[Trx + T h ] 


(4.6) 


and 


V. = G|T M + r«]. 


(4.7) 


The Yfactor is defined as 


v _Y±— + Th 

flctor - V e ~ Tnx + Tc' 


Solving for Trx with T k = 273 K and T e = 77 K yields 


(4.8) 


Trx — 


273-77 Y f>ctor 

^factor 1 


(4.9) 


The receiver temperatures (double side band) measured at 215 and 230 GHz were 
220 K and 240 K respectively. 

The primary objective of our observation was to detect the potential dip in 
Jupiter’s spectrum resulting from HjS absorption. Our approach was to use the 
total band width of the AOS essentially as a band pass filter or photometer at 
two or more frequencies. The emission is measured near the center of the HjS line 
where the potential absorption from HjS is at a maximum. The emission is then 
measured on either side (or both sides) of the line where the H 2 S opacity is small 
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the continuum ). In order to infer H,S absorption, we need only measure a 
differential emission between two or more wavelengths. Lellouch ct al. (1984) used 
a similar approach to search for PH S and HCN spectral features on Jupiter. A 
second objective of our observation was to obtain a reliable brightness temperature 
of Jupiter at this wavelength. This requires a precise estimate of the brightness 
temperature of Mars, which was used as the primary calibrator. 

Figure 4.3 shows the potential effect of HjS on Jupiter’s spectrum at 216 
GHz (1.4 mm). The dotted and solid lines are computed Jovian spectra with and 
without HjS absorption, respectively. We also show the frequencies at which we 
observed Jupiter with the double side band CSO receiver. We observed Jupiter 
with the local oscillator (LO) centered at 215.3 GHz and 229.6 GHz. We were 
unable to observe at the low frequency tail of the HjS line, because there were no 

available oscillators operating near 200 GHz. 

We observed both Jupiter and Mars on (UT) 25 November 1990 with the LO 
centered at 215.3 GHz. The frequency of the HjS line (216.7 GHz) is centered 
in the upper side band. We observed Mars first, then Jupiter, followed by Mars 
again, checking for variation in the observed temperature of Mars. We observed 

no variation in the temperature of Mars on this night. 

We used the same procedure to observe the two planets on (UT) 26 Novem- 
ber, 1990 at a frequency 229.6 GHz. This frequency was chosen so that the CO 
transition at 230.5 GHz, which has been observed in the spectrum of Mars, would 
be between the upper and lower side bands. Therefore, it would not interfere 
with the continuum of Mars. At this frequency, some variation occurred in the 
observed antenna temperature of Mars prior to and following the observation of 
Jupiter. This variation results in a large uncertainty in the observed temperature 
of Mars and thus a large uncertainty in the Mars/Jupiter ratio at this frequency. 
The receiver was then retuned to a center frequency of 215.3 GHz to repeat the 
observation of 25 November. We observed Mars followed by Jupiter. Because the 
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Figure 4.3: Effect of H]S absorption at 216 GHz and frequencies at which we 
observed Jupiter with the double side band (DSB) CSO receiver. Dashed line: 
spectrum with HjS absorption, solid line: spectrum without HjS absorption. 
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elevation of Mars was low and weather conditions marginal after observing Jupiter, 
we were unable to observe Mars a second tune. 

Ideally, the observations of Jupiter should be made when the calibrator 
(Mars) is close to the source (Jupiter) in the sky. In that case, the effects of 
any temporal or spatial variation in the earth’s atmospheric opacity would be lim- 
ited. However, good observations of Jupiter relative to Mars can still be made even 
if the planets are not ideally positioned. During our observation, Jupiter reached 
zenith approximately 5 hours after Mars. We were careful to observe Jupiter and 
Mars at times when the atmospheric opacity was observed to be stable and both 
planets were at similar elevations. The amount of observing time satisfying these 
conditions was limited. 

The airmass, M, is a term commonly used by astronomers to describe the 
elevation of the target. The airmass is related to the zenith angle 6 (90 degrees - 
elevation) by 

M = sec(0) = — . (4-10) 

H 

A plot of the airmasses of Jupiter and Mars is shown in Figure 4.4. Generally, ob- 
servations should be made when M < 1.5. Jupiter and Mars were both positioned 

at M < 1.5 between 02:00 and 04:00 local time. 

Chopping was accomplished by position switching or alternately pointing the 
telescope ON the source (planet) and then at a position in the sky OFF the source 
5' in either the + or - azimuth direction. The telescope remains both ON the 
source and OFF the source for a duration of 10 seconds. One scan of the source 
is defined as 4 ON /OFF cycles resulting in a total integration time of 80 seconds. 
The basic symmetric pattern is ON -OFF +OFF ON, where OFF is 5 off the 
planet in the + and - azimuth directions. Normally, the scans alternately begin in 
the ON and OFF positions. This alternating pattern may be manually overridden. 

The size of Jupiter is approximately equal to the beam width of the tele- 
scope. Therefore, telescope pointing is critical. Accurate telescope pointing is 


Airmass 
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accomplished by constructing a five point map of the planet. An on-line program 
uses a Gaussian fit to position the planet in the center of the telescope beam. 
For observations of Jupiter, the telescope was typically repositioned after every 
other scan. The telescope was repositioned following every third or fourth scan 
of Mars. The pointing remained accurate to better than 4 W (or ~ 10%) between 
repositioning. 


4.3 Calibration 


Before each scan of the planet, the computer automatically executes a chopper 
wheel calibration scan. This calibration procedure removes the effects of the 
earth’s atmospheric opacity. Because both our source and calibrator have finite 
and different angular sizes, additional calibration steps are required. 

During the chopper wheel calibration scan, the computer controls a chopper 
wheel which located in front of the Cassegrain feed (between the secondary and 
tertiary mirrors). First, a hot (ambient temperature) load is positioned in the 
beam of the antenna. The back-end computer inserts the necessary attenuators to 
prevent the AOS from saturating. A voltage level, V h , is measured with the AOS. 
The hot load is removed and the telescope is positioned on the sky (off source). A 
voltage, V^kyi is measured with the AOS. A calibration factor, Cal, is defined as 


_ Vi, - V.ky _ O « P 2 IhLl 

k»ky y*ky 


(4.11) 


(Phillips, 1989) where G is the system temperature to voltage conversion factor, 
a is called the hot spillover efficiency (1 - the fraction of power falling on the 
ground), /? is the cold spillover efficiency (1 - the fraction of power falling on the 
sky), 7 is the source coupling efficiency, and r is the atmospheric opacity at the 
zenith angle in optical depths. A similar ratio, OO-scan (On/ Off scan), is defined 
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for a scan of the source (planet) 


00 — scan = 




(4.12) 


V,„ V. kr 

where V, is the voltage level of the source. Taking the ratio of the two quantities 
OO — scam and Cal yields 

OO — scam PlT, , . 

- (4.13) 


or 


r,= 


Cal T h 

7\ OO — scan 


(4.14) 


01 Cal 

for double side-bamd observations. Since the vadue of /?7 is not known, the com- 
puter (CLASS data amalysis package) gives a plot of the corrected source temper- 
ature (TX) versus upper and lower side-bamd frequency where 

2 T* OO — scam 


T* — 

A M 


Cal 


(4.15) 


A typical Mars spectrum is shown in Figure 4.5. 

For a source (plamet) whose size is comparable to the size of the antenna beam, 
the product of the cold spillover efficiency, /?, amd the source coupling efficiency, 
7 , cam be computed from: 


/? 7 (MainBeam) = * [l ~ exp[- (^) In 2]] (4.16) 

(Phillips, 1989) where D is the angle subtended by the plamet, 6 is the beam 
width (FWHM), Jpianet is the actual brightness temperature of the plamet, amd T\ 
is the observed temperature (automatically corrected by the computer for atmo- 
spheric effects amd telescope losses). Mars is am ideal calibrator, because it is both 
bright amd compact. Because Mars has a tenuous atmosphere, its millimeter- wave 
spectrum is relatively flat (with the exception of a few narrow lines such as the 
CO line mentioned above). Its brightness temperature can be predicted with good 
accuracy. By assuming a value of T^ ar , based on radiative tramsfer models and 
using the observed value of TX (Mars), a value of 01 cam be computed. This value 
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108 


of ( S'! is then used to compute the brightness temperature of Jupiter, Tj up i ter , from 
the observed value T\ (Jupiter). 

In Equation 4.16, the second term represents a correction for the partial filling 
of the antenna beam (Ulich et ai ., 1980). We compute the size or diameter of the 
planets as viewed from the earth, 2?, (in arc seconds) from the precise earth to 
planet distance and the planet’s radius. For Jupiter, we use an equatorial radius 
at the 1 bar level of 71495 km and an ellipticity of 0.065. This value is based on 
the radio occultation experiment aboard Voyager (Lindal et a!., 1986). We use a 
value of 3397 km for the Martian equatorial radius with an ellipticity of 0.0006. 
These are the same values as those used by Griffin et ai. (1986). The exact earth to 
planet distances for Mars and Jupiter were taken from the Astronomical Almanac 
(1990). For Mars and Jupiter, the true geocentric distances were approximately 
4.86 and 0.4 AU respectively. This yields a diameter of ab out 18” for Mars and 
polar and equatorial diameters of 37.9" and 40.5" respectively for Jupiter. The 
effective diameter of Jupiter in Equation 4.16 is the geometric mean of the polar 
and equatorial diameters. 

The beam width of the antenna is inversely proportional to the wavelength 


A. The FWHM beam width of a parabolic reflector is 

(417) 

(Gagliardi, 1984) where d is the diameter of the reflector, and p a is the antenna 
efficiency factor which is dependent on how the antenna aperture is electromag- 
netically illuminated. Because p a is not precisely known, the beam width should 
be measured. In order to infer the beam width of the antenna (0), we obtained 
a map of Jupiter at 230 GHz from Schinckel ( private communication , 1991). The 
contour plot is shown in Figure 4.6. We read off values from the two dimensional 
beam map of Jupiter at the equator (using linear interpolation between the con- 
tour lines) to develop a one dimensional beam map of Jupiter’s equator. The map 
of Jupiter’s equator was then deconvolved with a model of the expected emission 
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at Jupiter’s equator which takes into account the effects of its oblate shape and 
limb darkening (see Section 3.1). Figure 4.7 shows the equatorial beam scan and 
the expected brightness temperatures of Jupiter’s equator. The resulting antenna 
pattern is shown in Figure 4.8. The beam is Gaussian with a FWHM of 45.9" at 
230 GHz. Assuming that the antenna or illumination efficiency factor is constant 
over the range of 215.3 to 230 GHz, the FWHM at 215.3 GHz would be 49.0". 

We use brightness temperatures of 211.9 and 212.6 K for Mars at 216 and 
230 GHz, respectively (Rudy et o/., 1987 and Rudy, private communication, 1991). 
Rudy’s model accounts for both seasonal and longitudinal variation. The longitu- 
dinal variation in the Martian brightness temperature predicted by Rudy’s model 
was less than 1 K or 0.5%. Since Rudy’s model is extrapolated from observations 
at 2 and 6 cm, we have used an additional check. Ulich (1981) has developed 
the only empirical model for predicting the millimeter brightness temperature of 
Mars. Ulich’s model is based on millimeter observations using absolute calibra- 
tion. Ulich’s model is limited in that it is based on observations at a single fre- 
quency (86 GHz or 3.5 mm). His model takes into account only seasonal variation 
and neglects longitudinal variation. At this frequency Rudy’s model predicts the 
brightness temperatures over one Mars rotation to vary between 197 and 203 K. 
Ulich’s empirical expression is 

Tb (86 GHz)=n(-^) / , ( 4 - 18 ) 

where R is the heliocentric distance, T B = 206.8±1.5 K, and R„ = 1.524 AU. 
During our observation, R = 1.50545 AU. The resulting Martian brightness tem- 
perature using Equation 4.18 is 208.1 ± 1.5 K at 86 GHz. This temperature is 
withiu about 5% of Rudy’s predicted brightness temperatures. Previous observers 
have arbitrarily assumed the uncertainty in the Martian brightness temperature 
to be ±10% or approximately ±20 K (see, e.g., Griffin et a/., 1986). Although 
Rudy estimates that the uncertainty in his model is better than 10% (Rudy, pri- 
vate communication, 1991), we will conservatively assume the uncertainty in the 
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Figure 4.7: Equatorial cut through the beam map of Jupiter and the computed 
brightness temperature (normalized to 1) for an equatorial cut through Jupiter 
including limb darkening 
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Figure 4.8: Plot of the beam shape of the CSO at 230 GHz; Solid line: Inferred 
beam shape of the CSO, dashed line: Gaussian beam with FWHM = 45.9*. 
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Martian brightness temperature to be 10%. 


4.4 Atmospheric Conditions 

The atmospheric opacity was monitored with a separate tipping radiometer or tau 
meter operated by the NRAO at 225 GHz. Plots of r. (the zenith optical depth) 
for 25 and 26 November are shown in Figures 4.9 and 4.10, respectively. These 
plots are good indicators of the atmospheric opacity and variability of the weather. 

The tau meter records opacity as a function of altitude in only one direction. 
We checked the weather conditions frequently for local variation (».«•, clouds). 
A 6 inch optical telescope (Paigeatron) mounted just below the CSO reflector 
and connected to a video display in the control room allowed us to monitor the 
visibility of the planets. Frequent trips outside the observatory were also required 
to observe local conditions. 

Another indicator of the observing conditions is the system temperature [T , yi ) 
which is defined as the system noise at the source elevation divided by the loss 

factor of the source or 

_ Trx + (1 ~ ae r ) r *> (4.19) 

lfyt ~ afiie-* 

(Phillips, 1989). The voltages of the source (planet) and the sky (V, and V. ky , 
respectively) are related to the receiver temperature (Trx) by 

V, ky = G[Trx + (1 - oce'^Tk] (4-20) 

and 

V. = G{Tnx + (1 - ae- T ^)T H + aW^), (4.21) 

where T h is the temperature of the sky, and T, is the temperature of the source. 
Substituting Equation 4.20 into Equation 4.11 yields 

ae- r I\ 

Trx + (1 — ae -r )T\ 


Cal = 


(4.22) 
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Then, T^, can be written as 


T _ T h 
I ' y *~ Cal^' 


(4.23) 


Because the quantity /?7 is not known, the computer gives the system temperature 
as 

^ = CaP ( 4 * 24 ) 

The system temperature for each scan is listed in Tables 4. 1-4 .3. 

We were unable to observe on the first scheduled night of observation (24 
November) due to poor weather conditions (low lying cumulus cloud and fog). We 
utilized this night by familiarizing ourselves with the receiver tuning process. The 
zenith optical depth (r # ) was stable at 0.05 on 25 November for approximately 
three hours. As the night progressed, the humidity rose from an initial 20% to 
100% at around 04:00 local time. At this time, fog was visible and observations 
were terminated. On 26 November, t x was higher but relatively stable at 0.2 for 
approximately four hours. The humidity was initially high (80%) and rose to 100% 
at 05:00. Fog was again visible, and we noted ice on the surface of the telescope 
at this time. 


4.5 Data Analysis 

Using the data analysis package CLASS, we fit a zero order baseline to the observed 
spectrum of each scan. The antenna temperature (TJ) corrected for atmospheric 
opacity is averaged over the 1024 channels of the AOS to reflect a mean and RMS 
deviation for in each scan. Tables 4.1-4.3 list T\, the RMS noise, the system 
temperature (T^,), and the zenith optical depth (r) for each scan. Also listed is 
the position of the antenna at the beginning of the scan (ON or OFF source). We 
also note when the antenna was repositioned. 

The overall observational uncertainty was dominated by scatter in the means 
of individual scans and not the base line ripple observed in a single scan. As 
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Table 4.1: Listing of observations (UT) 25 November 1990 at 216 GHz. 

Source T\ (K) RMS (K) T„. r OW/ OFF Comment! 

Not pointed 


47.5 3.1 

43.2 1.5 


1014 

1016 

Cel 

Men 

50.4 

1.5 

063.2 

0.055 

OFF 

1016 

1017 

Cel 

Men 

58.4 

1.4 

050.2 

0.056 

ON 

1018 

1010 

Cel 

Men 

50.4 

1.3 

963.2 

0.055 

OFF 

1020 

1021 

Cel 

Men 

58.7 

1.3 

070.2 

0.055 

ON 

1022 

1023 

Cel 

Men 

50.3 

1.4 

066.7 

0.055 

OFF 

1024 

1025 

Cel 

Men 

60.4 

1.2 

080.6 

0.055 

ON 

1026 

1027 

Cel 

Men 

59.5 

1.7 

076.0 

0.055 

OFF 

1028 

1020 

Cel 

Men 

60.3 

1.6 

984.0 

0.056 

OFF 

1030 

1031 

1032 

1033 

1034 

1035 

Cel 

Cel 

Cel 

Cel 

Cel 

Jupiter 

162.8 

2.7 

1013 

0.055 

ON 

1036 

1037 

Cel 

Jupiter 

148.2 

3.5 

1013 

0.055 

OFF 

1038 

1030 

Cel 

Jupiter 

150.1 

1.0 

091.3 

0.055 

ON 

1040 

1041 

Cel 

Jupiter 

149.9 

2.0 

096.5 

0.055 

OFF 

1042 

1043 

Cel 

Jupiter 

162* 

2.2 

096.1 

0.056 

OFF 

1044 

1045 

Cel 

Jupiter 

146.7 

3.1 

086.1 

0.055 

ON 

1041 

1048 

Cel 

Men 

57.8 

0.8 

062.5 

0.055 

ON 

1040 

1050 

Cel 

Men 

62.0 

0.0 

077.8 

0.055 

ON 

1051 

1052 

Cel 

Men 

61.1 

0.0 

074.3 

0.055 

OFF 

1053 

1054 

Cel 

Men 

61.7 

1.1 

003.0 

0.056 

OFF 

1055 

1056 

Cel 

Men 

57.5 

0.0 

080.3 

0.055 

OFF 


2 Five-point* 


1 Five-point 


3 Five-point* 


1 Five-point 


1 Five-point 


1 Five-point 


1 Five-point 
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Table 4.2: Listing of observations (UT) 26 November 1990 at 230 GHz. 


Sean Source TjJ (K) RMS (K) T. t . r ON/OFF Comment* 
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Table 4.3: Listing of observations (UT) 26 November 1990 at 216 GHz. 


Scan 

Source 

n ( k ) 

RMS (K) 

Tgyg 

r 

ON/OFF 

Comments 

1057 

Cal 






1 five point 

1058 

Mars 

47.5 

0.8 

1512 

0.226 

ON 


1059 

Cal 







1060 

Mars 

43.2 

0.7 

1540 

0.224 

OFF 


1061 

Cal 







1062 

Mars 

40.8 

0.9 

1559 

0.224 

ON 


1063 

Cal 







1064 

Mars 

49.2 

1.5 

1591 

0.224 

OFF 


1065 

Cal 






2 Five-points 

1066 

Jupiter 

129.9 

2.2 

1093 

0.224 

OFF 


1067 

Cal 







1068 

Jupiter 

144.0 

2.8 

1113 

0.230 

ON 


1069 

Cal 






2 Five-point 

1070 

Jupiter 

132.5 

2.2 

1113 

0.236 

OFF 


1071 

Cal 






2 Five-points 

1072 

Jupiter 

147.0 

2.4 

1117 

0.236 

ON 


1073 

Cal 






2 Five-point 

1075 

Cal 







1076 

Cal 







1077 

Orion 







1078 

Cal 







1079 

Orion 







1080 

Cal 







1081 

Orion 







1082 

Cal 







1083 

Orion 







1084 

Cal 






1 Five-point 

1085 

Jupiter 

147.7 

1.7 

1079 

0.221 

ON 


1086 

Cal 







1087 

Jupiter 

137.3 

3.4 

1083 

0.218 

OFF 


1088 

Cal 






1 Five-point 

1089 

Jupiter 

147.4 

1.9 

1077 

0.218 

ON 


1090 

Cal 







1091 

Jupiter 

139.3 

3.1 

1082 

0.219 

OFF 


1092 

Cal 






1 Five-point 

1093 

Jupiter 

151.5 

1.7 

1085 

0.219 

ON 
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expected, we are limited by 8ky noise and systematic errors, not the sensitivity 
of the receiver. The averaged observed temperatures, r;, are listed in Table 4.4. 
We also list the computed values of /?7. The expected value of ^7 at 230 GHz is 
approximately 0.72. We note that our observed values of /?7 varied from night to 
night by as much as 10%. 

The spectra of both Jupiter and Mars at the two observed frequencies are 
expected to be flat over band width (500 MHz) of the receiver. Both random and 
systematic processes contribute to the observed baseline ripple (see Figure 4.8). 
The ripple is still present after averaging several scans. Therefore, it must be 
systematic in part. 

We also note systematic variations in the scans of both Jupiter and Mars. 
The recorded temperatures of Jupiter from scans which began with the telescope 
ON the planet were consistently higher than those from scans which began with the 
telescope OFF the planet (on the sky). We believe that the computer-controlled 
attenuators were not set to adequate levels on the scans which began OFF the 
source. This would cause the amplifiers to be driven into compression which would 
result in lower observed antenna temperatures. Compression of the amplifiers was 
also observed and confirmed independently by Schinckel ( private c ommun ication, 
1991). Therefore, we believe that only the scans which started ON Jupiter are reli- 
able and include only these scans in our calculation of Jupiter’s observed brightness 
temperatures. We list the data used to derive brightness temperatures for Jupiter 
at 215 and 230 GHz in Table 4.5, The uncertainty in the relative brightness tem- 
peratures of Jupiter, T Jupit „, for individual nights is the quadratic sum of the 
observed uncertainties in T\ for Mars and Jupiter (note: This does not include 
the additional ±10% uncertainty in the absolute temperature of Mars). 

Our observed brightness temperatures of Jupiter and those of previous ob- 
servations near 1.4 mm are reported in Table 4.6. The first uncertainty represents 
the 1 o statistical uncertainty in the Mars/Jupiter ratio appropriate for searching 
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Table 4.4: Observed antenna temperatures. 


Freq. Date Source T\ ± ler 07 Comments 

(GHz) (UT) (K) 


215.31 

11/25/90 

Mars 

59.7±0.4 

215.31 

11/25/90 

Jupiter 

154.9±6.7 

161.6±1.2 

148.2±0.9 

215.31 

11/26/90 

Mars 

51.4±0.6 

215.31 

11/26/90 

Jupiter 

141.8±6.4 

147.5±1.2 

134.7±2.5 

229.60 

11/26/90 

Mars 

60.5±4.4 

54.8±0.4 

63.6±0.7 

229.60 

11/26/90 

Jupiter 

154.5±12.9 

169.2±1.5 

143.4±1.3 


0.825 All scans (2 sets) 

All scans (1 set) 
scans starting ON source only 
wfjwiw starting OFF source only 
0.713 All scans (1 set) 

All scans (2 sets) 
scans starting ON source only 
scans starting OFF source only 
0.743 All scans (2 sets) 

0.673 Low’s only 

0.782 High’s only 

All scans (1 set) 

scans starting ON source only 

scans starting OFF source only 
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Table 4.5: Observational data used to computed Jupiter’s brightness temperature. 


Freq. 

(GHz) 

Date 

(UT) 

TX±la 

Mars 

Tl±lc 

Jupiter 

assumed 

T. J u piter 
observed 

215.31 

11/25/90 

59.7±0.4 

161.6±1.2 

211.9 

166.3 ±1.7 

215.31 

11/26/90 

51.4±0.6 

147.5±1.2 

211.9 

175.0 ± 2.5 

229.60 

11/26/90 

60.5±4.4 

169.2±1.5 

212.6 

178.1 ± 13.0 

229.60 

11/26/90 

54.8±0.4 

169.2±1.5 

212.6 

196.7 ± 2.3 

229.60 

11/26/90 

63.6±0.7 

169.2±1.5 

212.6 

169.4 ± 2.4 


for possible spectral features. The second uncertainty includes an additional 10% 
uncertainty in the assumed brightness temperature of Mars required for absolute 
brightness temperature measurements. Our observed brightness temperatures are 
in good agreement with previous observations at wavelengths near 1.4 mm (Rather 
et al.y 1974, Courtin et al ., 1977, Ulich et al ., 1984 and Griffin et al.> 1986). All 
of the previous observations were made with broadband filters (see also Table 3.3) 
and only the observations of Courtin et al. (1977) and Griffin et al. (1985) used 
Mars as a calibrator. Rather et al. (1974) and Ulich et al. (1984) somewhat arbi- 
trarily assigned brightness temperatures of 150 and 165 K, respectively, to Jupiter 
based on their observations of several planets and previously observed planetary 
brightness temperatures at longer and shorter wavelengths. 

Finally, we chose to utilize the spectral information available in the 1024 
channels of the AOS to obtain additional calibration information and to check for 
the unlikely presence of narrow line core emission from HjS in Jupiter’s strato- 
sphere. We observed line emission from the core of the Orion molecular cloud on 
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Table 4.6: Observed Jovian brightness temperatures near 1.4 mm. 


Wavelength 

(mm) 

Freq 

(GHz) 

T| 

(K) 

1.30 

231 

175.0±2.5(18) 

1.32 

227 

170.9±3.9(18) 

1.32 

227 

165±8(18) 

1.39 

216 

178.1 1±13.0(22) 

1.40 

214 

148±16(22)‘ 

1.40 

214 

168±ll(20) e 


Calibrator Ai/ Reference 

(GHz) 


Mars 

3.2 

This work 

Mars 

70 

Griffin et o/.(1986) 

planets 

39 

Ulich et a/. (1984) 

Mars 

3.2 

This work 

planets 

275 

Rather et a/. (1975) 

Mars 

210 

Courtin et a/.(1977) 


fl First uncertainty is l<r observational uncertainty, 

the second uncertainty includes an additional 10% for absolute calibration 
‘Brightness temperatures as corrected in Berge and Gulkis (1976) 

* Recalculated using beam correction factor in Ulich (1980) 
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both nights at 215.3 GHz in addition to the two planets. The hydrogen sulfide line 
emission from Orion at 216.7 GHz (centered in the upper side band) was clearly 
visible after one scan. If an HjS emission core is present in the scans of Jupiter, 
it will also appear in the center of the upper side band. An H 2 S mixing ratio of 
close to 1 x 10 -6 would have to be present in the stratosphere in order to detect 
an emission core. Based on the strict upper limits (as low as 1 x 10"®) which 
have been placed on the abundance of HjS in Jupiter’s stratosphere at several 
wavelengths from the infrared to uv (summarized by Larson et al., 1984), no such 
feature is expected. Figure 4.11 shows an average of all scans of Jupiter at 215.3 
GHz for 25 November. Indeed, no H 2 S emission core is visible. 

The line intensities in the spectra of Orion also provide a secondary calibra- 
tion source. An averaged spectrum of Orion observed on 26 November is shown in 
Figure 4.12. There is good agreement between our spectra of Orion (26 November) 
said those of Sutton et al. (1985). However, the line intensities of the averaged 
spectra taken on 25 November are significantly higher than those of Sutton et al. 
(1985) which could be due to a side band ratio ^ 1. This might also explain the 
high value of /?7 observed on 25 November. Therefore, we believe the 215.3 GHz 
observation on 26 November to be more reliable than the 25 November observation. 


4.6 Observational Results 

Figure 4.13 shows several synthetic emission spectra of Jupiter near 216 GHz using 
the radiative transfer and thermochemical model described in Chapter 3. We have 
computed the Jupiter’s emission using the NH S and H 2 S vertical distributions 
shown in Figure 3.9. Our observed brightness temperatures and corresponding 
error bars are also shown along with several other previous observations. 

We were unable to obtain repeatable observations at individual frequencies 




Figure 4.11: Average of Jupiter scans at 215.3 GHz for 25 


November 


213700 213800 213900 214000 214100 

Image Frequency (MHz) 



®07S$-NeDU-aiB90H H2S22 0_2 1 1 CSO 500 MHz O: 26-NOU-199 

: 0 5:32:03.000 DEC: -5:24:21.00 (1950.0) CH * s : 0.0 

Ertown23 . 73u : 0.2210 T sys : 1693. Time: 3.235 
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Figure 4.12: Average of Orion scans at 215.3 GHz for 26 November 
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Jupiter's observed and computed spectrum near 1.4 mm 



Figure 4.13: Observed spectrum of Jupiter near 216 GHz. •: our observations; 
o: other observations from Table 4.6. Theoretical spectra using distributions in 
Figure 3.9, solid line: NH 3 , HjO, and pressure-induced opacity only; dashed line: 
HjS opacity added. 
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over several nights due to poor weather conditions. The uncertainties for our ob- 
servation (especially for the single observation at 229.6 GHz) are large. Therefore, 
we did not detect a statistically significant difference in emission due to H 2 S at 
the two wavelengths. Although we were unable to set tighter limits on the H 2 S 
abundance in Jupiter’s atmosphere, our results reinforce previous limits summa- 
rized by Larson tt al., 1985. Even though we did not achieve our primary objective 
of detecting H 2 S gas, we were able to measure Jupiter’s brightness temperature 
relative to Mars at the highest spectral resolution recorded at this wavelength. 
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CHAPTER 5 

Summary and Conclusion 

5.1 Uniqueness of Work 

This research makes contributions in the areas of both planetary science and 
millimeter-wave measurements. The major contributions are su mmar ized below: 

1. Conducted first laboratory measurements of gaseous NH S opacity at Ka-band 
(u = 32- 40 GHz, A = 7.5 — 9.38 mm) under simulated Jovian conditions 

2. Conducted first laboratory measurements of gaseous NH S opacity at W-band 
[y = 94 GHz A = 3.2 mm) under simulated Jovian conditions 

3. Developed new parameterization of the Ben-Reuven line shape for computing 
the absorptivity of NHs under Jovian conditions 

4. Conducted first laboratory measurement of the hydrogen and helium broad- 
ened line width of gaseous HjS at 1.4 mm 

5. Documented for the first time new techniques in millimeter-wave measure- 
ments (t.g., laser tuning of the resonator) and documented for the first time 
solutions to problems which occur in millimeter-wave measurements (c.g., 
dielectric loading) 

6. Incorporated for the first time the results of the millimeter-wave laboratory 
absorption measurements of NHj and HjS into a radiative transfer model 
and computed synthetic emission spectra from the four giant planets 
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7. Developed new expressions for computing the radio opacity of pressure- 
induced absorption from Hj-Hj, Hj-He and Hj-CH^ pairs and water vapor 
under Jovian conditions 

8. Demonstrated that Jupiter’s observed millimeter-wave emission is entirely 
consistent with synthetic Bpectra in which NHj is the primary source of 
opacity 

9. Observed Jupiter at 1.4 mm for the first time using a high resolution receiver 
and obtained a reliable brightness temperature of Jupiter using Mars as the 
calibration standard 

This research has culminated in the publication of several journal papers: 
Joiner et al. (1989), Joiner and Steffes (1991), and Joiner et al. (1991). This 
research has also been presented as several conference reports: Joiner et al. (1987), 
Joiner and Steffes (1988), Joiner and Steffes (1989), Joiner and Steffes (1990a, b). 

5.2 Suggestions for Future Research 

Although we were unable to obtain a positive detection or HjS or set new up- 
per limits on the H 2 S abundance in Jupiter’s atmosphere, the search for H 2 S at 
millimeter wavelengths should not be abandoned. Poor weather conditions and in- 
strumental difficulties contributed to large uncertainties in our observation. How- 
ever, given more ideal conditions, the statistical uncertainties could be significantly 
reduced and systematic effects eliminated. In order to reduce the observational 
uncertainty, the integration time must be increased (i.e., observations should be 
repeated over several days). Weather conditions must be exceptional in order to 
eliminate any variations in the observed temperatures of Jupiter and Mars caused 
by variations in the earth’s atmospheric opacity. In addition, observations should 
be made at a minimum of three frequencies, corresponding to the center of the 
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HjS line and on both sides of the line, in order to achieve a convincing detection 
of the pressure-broadened line. 

Improvements in millimeter- wave technology will make the detection of broad 
spectral features (i. e., pressure-broadened lines in planetary atmospheres) more 
feasible in the future. Wide band oscillators coupled with computer-controlled 
tuning will significantly reduce the time required to tune receivers. In addition, 
new instruments which can be mounted on existing millimeter and submillimeter 
antennas, such as Fourier Transform Spectrometers, will be better suited for the 
detection of broad spectral features. In the future, spectrometers with the proper 
resolution (on the order of GHz) for planetary spectroscopy will be utilized to 
improve upon the results obtained with the multi-wavelength approach. Improved 
calibration techniques will reduce the large uncertainty that presently exists in the 
absolute calibration of millimeter planetary observations. 

Although this research represents great strides in the understanding of the 
millimeter- wave spectra of the Jovian planets, more microwave and millimeter- 
wave laboratory measurements are needed in order to fully interpret all of the 
available planetary observations. For example, measurements of the millimeter 
dielectric properties of solid NHj and NH4SH are needed to accurately assess 
the potential effect of cloud condensates on the millimeter-wave spectra of the 
giant planets. Absorption measurements of pressure-broadened H 2 S at centimeter 
wavelengths are needed in order to interpret the centimeter emission from Uranus 
and Neptune. In addition, a measurement of the equilibrium constant for the 
reaction of NH S and H 2 S to form solid NH4SH is needed in order to accurately 
predict the altitude at which the putative NH4SH cloud will form. It is hoped that 
this research has provided the initiative for continued exploration of the millimeter- 
wave spectrum as part of the study of the Jovian planets. 




132 


BIBLIOGRAPHY 


ASTRONOMICAL ALMANAC 1990. U. S. Naval Observatory, U. S. Govt. Print- 
ing Office, Washington, D. C. 

ANDERS, E. AND N. GREVESSE 1989. Abundances of the elements: mete- 
oritic and solar. Geochim. Cosmochxm . Acta 53, 197-214. 

ATREYA, S. K.AND P. N. ROMANI 1985. Photochemistry and clouds of 
Jupiter, Saturn, and Uranus. In Recent advances in planetary meteorology. 
ed. G. E. Hunt. Cambridge Univ. Press, Cambridge, 17-68. 

BAARS, J. W. M., R. GENZEL, I. I. PAULINY-TOTH, AND A. WITZEL 
1977. The absolute spectrum of Cas A: an accurate flux density scale and 
a set of secondary calibrators. Astron. Astrophys. 61 , 99-106. 

Bachet, G., E. R. Cohen, P. Dore, and G. Birnbaum 1983. The 

translational-rotational absorption spectrum of hydrogen. Can. J. Phys 61 , 
591-603. 

BATTAN, L. J. 1973. Radar Observation of the Atmosphere. University of 
Chicago Press, Chicago. 

BEN-REUVEN, A. 1966. Impact broadening of microwave spectra. Phys. Rev. 
145, 7-22. 

BERGE, G. L., AND S. GULKIS 1976. Earth-based radio observation of 
Jupiter: Millimeter to meter wavelengths. In Jupiter (T. Gehrels, Ed.), 621- 
692. Univ. of Arizona Press, Tuscon. 

BEZARD, B., A. MARTEN, J. P. BALUTEAU, D. GAUTIER, J. M. FLAUD, 
AND C. CAMY-PEYRET 1983. On the detectability of H 2 S in Jupiter. 
Icarus 55, 259-271. 

BJORAKER, G. L., H. P. LARSON, AND V. G. KUNDE 1986A. The gas 
composition of Jupiter derived from 5 pm airborne spectroscopic observa- 
tions. Icarus 66, 579-609. 

BJORAKER, G. L., H. P. LARSON, AND V. G. KUNDE 1986B. The abun- 
dance and distribution of water vapor in Jupiter’s atmosphere Astrophys. J. 
311 , 1058-1072. 



133 


BORYSOW, J., L. TRAFTON, AND L. FROMMHOLD 1985. Modeling of 
Pressure- Induced Fax-Infrared Absorption Spectra: Molecular Hydrogen Pairs. 
Astrophys. J. 296, 644-654. 

BRIGGS, F. H., AND P. D. SACKETT 1989. Radio observations of Saturn 
as a probe of its atmosphere and cloud structure. Icarus 80, 77-103. 

CRC PRESS 1983. CRC Handbook of Chemistry and Physics (R. C. Weast, 
Ed.). CRC Press, Boca Raton, FL. 

Chedin, A., N. Husson, N. A. Scott, I. Cohen Hellalet, A. Berroir 
1982. GEISA data bank (No 3). 

Cogdell, J. R., J. H. Davis, B. T. Ulrich, and B. J. Wills 1975. 
Flux density measurements of radio sources at 2.14 millimeter wavelength. 
Astrophys. J. 196, 363-368. 

COLLIN, R. E. 1966. Foundations for Microwave Engineering. McGraw-Hill, 
New York. 

COURTIN, R., N. CORON, T. ENCRENAZ, R. GlSPERT, P. BRUSTON, 

J. LEBLANC, G. DAMBIER, A. VIDAL- MADJAR 1977. Observations 
of giant planets at 1.4 mm and consequences on the effective temperatures. 
Astron. Astrophys. 60, 115-123. 

DEPATER, I., AND S. T. MASSIE 1985. Models of the millimeter-centimeter 
spectra of the giant planets. Icarus 62, 143-171. 

DEPATER, I. 1986. Jupiter’s Zone Belt Structure at Radio Wavelengths: n. 
Comparison of Observations with Atmosphere Calculations. Icarus 68, 344- 
365. 

DEPATER, I., P. N. ROMANI AND S. K. ATREYA 1989. Uranus deep 
atmosphere revealed. Icarus 82, 288-313. 

DEPATER, I., P. N. ROMANI AND S. K. ATREYA 1991. Possible microwave 
absorption by HjS gas in Uranus’ and Neptune’s atmospheres. Icarus 91, 
220-233. 

DEPATER, I. AND M. RICHMOND 1989. Neptune’s microwave spectrum from 
1 mm to 20 cm. Icarus 80, 1-13. 

DEIRMENDJIAN, D. 1969. Electromagnetic Scattering on Spherical Polydis- 
persions , American Elsevier Publishing Co, Inc., New York. 



134 


DORE, P., L. NENCINI, AND G. BlRNBAUM 1983. Far infrared absorption 
in normal Hjfrom 77 to 298 K. J. Quant . Spec. Rad. Trans 30, 245-253. 

FLAUD, J. M., C. CAMY-PEYRET, AND J. W. C. JOHNS 1983. The far 
infrared spectrum of hydrogen sulfide. Can. J. Phys. 61, 1462-1473. 

GAGLIARDI, R. M. 1984. Satellite Communications , Lifetime Learning Pub- 
lications, Wadsworth, Inc., Belmont, CA. 

GASIEWSKI, A. J. 1990. Microwave radiative transfer in hydrometers. In 
Atmospheric Remote Sensing by Microwave Radiometry (M. A. Janssen, ed.) 
To be published by John Wiley and sons. 

GAUT N. E. AND E. C. REIFENSTEIN III 1971. Environmental Research 
and Technology, Inc., Report No. 13, Lexington, MA. 

GOODMAN, G. C. 1969. Models of Jupiter’s atmosphere Ph. D. thesis , Uni- 
versity of Illinois, Urbana. 

Griffin, M. J., P. A. R. Ade, G. S. Orton, E. I. Robson, W. K. 
Gear, I. G. NOLT, AND J. V. RADOSTITZ 1986. Submillimeter and 
millimeter observations of Jupiter. Icarus 65, 244-256. 

GROSS, E. P. 1955. Shape of collision-broadened spectral lines. Phys. Rev. 
97, 395-403. 

GROSSMAN, A. W. 1990. Microwave imaging of Saturn’s deep atmosphere 
and rings. Ph. D. thesis , California Institute of Technology, Pasadena. 

GULKIS, S. 1987. Radio astronomy, planetary. Encycl. Phys. Sci. Tech.,N 61. 
n, 633-654. 

GULKIS, S. AND I. DEPATER 1984. A review of the millimeter and centimeter 
observations of UranuB. In Uranus and Neptune : Proc. of a workshop held 
in Pasadena, CA, Feb. 6-8, 1984 , NASA Conf. Publ. 2330, 225-262. 

HELMINGER, P. AND F. C. DELUCIA 1972. Pressure broadening of hydrogen 
sulfide. J. Quant. Spectrosc. Rad. Trans. 17,751-754. 

HOFSTADTER, M. D., G. L. BERGE AND D. O. MUHLEMAN 1990. 
Vertical motions in the Uranian atmosphere: an analysis of radio observation. 
Icarus 84, 261-267. 



135 


HOFSTADTER, M. D. AND D. O. MUHLEMAN 1988. Latitudinal varia- 
tions of ammonia in the atmosphere of Uranus: an analysis of microwave 
observations. Icarus 81, 396-412. 

JOINER, J., JENKINS, J. M., AND STEFFES, P. S., (1987). Laboratory 
Measurements of the Opacity of Gaseous Ammonia (NH3) in the 7.3 - 8.3 tutti 
(36-41 GHz) Range Under Simulated Conditions for the Jovian Atmospheres, 
Bull , Amcr. Astron. Soc. 19, 838. 

JOINER, J., Jenkins, J. M., AND STEFFES, P. S., (1988). Millimeter-wave 
Measurements of the Opacity of Gaseous Ammonia (NHj) Under Simulated 
Conditions for the Jovian Atmospheres, Bull. Amer. Astron. Soc. 20, 867. 

JOINER, J., P. G. STEFFES, AND J. M. JENKINS 1989. Laboratory 
measurements of the 7.5-9.38 mm absorption of gaseous ammonia (NHj) 
under simulated Jovian conditions. Icarus 81, 386-395. 

JOINER, J., AND STEFFES, P. G. (1989). Models of the Millimeter-wave 
Emission of the Jovian Atmosphere Utilizing Laboratory Measurements of 
Gaseous Ammonia (NH S ). Bull. Amcr. Astron. Soc. 21, 945. 

JOINER, J., AND STEFFES, P. G. (1990). Study of the Millimeter- wave Ab- 
sorbing Constituents in the Jovian Atmospheres. Bull. Amcr. Astron. Soc. 
22, 1032. 

JOINER, J., AND STEFFES, P. G. (1990). Measurement of the Pressure- 
Broadened Linewidths of Gaseous Hydrogen Sulfide (HjS). Poster presented 
at Second International Conference on Laboratory Research for Planetary 
Atmospheres. 

JOINER, J., AND P. G. STEFFES 1990. Modeling of Jupiter’s millimeter- 
wave emission utilizing laboratory measurements of ammonia (NH S ) opacity. 
JGR:Plancts , accepted. 

JOINER, J., P. G. STEFFES, AND K. S. Noll 1990. Search for sulfur 
(HjS) on Jupiter at millimeter wavelengths. Submitted to IEEE Trans. Mi- 
crowave Thcor. Tech. 


Klein, M. j. AND S. GULKIS 1978. Jupiter’s atmosphere: Observations and 
interpretation of the microwave spectrum near 1.25 cm wavelength. Icarus 
25,44-60. 



136 


KLEIN, M. J., M. A. JANSSEN, S. GULKIS AND E. T. OLSEN 1978. 
Saturn ’8 microwave spectrum: Implications for the atmosphere and the rings. 
In The Saturn System , NASA CF-2068, (D. M. Hunten and D. Morrison, 
Eds.), 195-216. NASA, Washington, D. C. 

KRAUS, G. F., J. E. ALLEN, JR., AND L. C. COOK 1989. Vapor 
pressure measurements of hydrogen sulfide. Bull. Amer. Astron. Soc. 21 
948. 


LARSON, H. P., D. S. DAVIS, R. HOFMAN, AND G. L. BJORAKER 
1984. The Jovian atmospheric window at 2.7 microns: A search for HjS. 
Icarus 60, 621-639. 

LELLOUCH, E., T. ENCRENAZ AND F. COMBES 1984A. The observability of 
HCN on Jupiter in the millimeter range. Astron. Astrophys. 135, 365-370. 

LELLOUCH, T. ENCRENAZ AND E., F. COMBES AND P. DROSSART 1984B. 
The observability of HCN and PHj on Saturn in the millimeter range. As- 
tron. Astrophys. 135, 371-376. 

LELLOUCH, E., F. COMBES, AND T. ENCRENAZ 1984c. Microwave observa- 
tions of Jupiter and Saturn. Astron. Astrophys. 140, 216-219. 

LELLOUCH, E., J. L. DESTOMBES 1985. Search for minor atmospheric 
species in the millimeter range of Jupiter and Saturn. Astron. Astrophys. 
152, 405-412. 

LINDAL, G. F., G. E. WOOD, G. S. LEVY, J. O. ANDERSON, D. N. 
SWEETNAM, H. B. HOTZ, B. J. BUCKLES, D. P. HOLMES, P. E. 
DAMS, V. R. ESHLEMAN, G. L. TYLER, AND T. A. CROFT 1981. 
The atmosphere of Jupiter: An analysis of the Voyager Radio Occultation 
measurements. J. Geophys. Res. 86, 8721-8727. 

LINDAL, G. F., D. N. SWEETNAM, AND V. R. ESHLEMAN 1985. The 
atmosphere of Saturn: An analysis of the Voyager Radio Occultation mea- 
surements. J. Geophys. Res. 90, 1136-1146. 

LINDAL, G. F., J. R. Lyons, D. N. Sweetnam, V. R. Eshleman, 
D. P. HINSON, AND G. L. TYLER 1987. The atmosphere of Uranus An 
analysis of the Voyager Radio Occultation measurements. J. Geophys. Res. 
92, 14987-15002. 

MlE, G. 1908. Beitrage zur optik truber medien, speziell kolloidaler metala- 
sungen. Ann. Physik t 25, 377. 



137 


MORRIS, E. C., AND R. W. Parsons 1970. Microwave absorption by gas 
mixtures at pressures up to several hundred bars. Astron. J. Phys. 23, 
335-349. 

PHILLIPS, T. G. 1989 Caltech Submillimeter Observatory (CSO) Handbook. 

PALTRIDGE, G. W. AND C. M. R. PLATT 1976. Radiative Processes in 
Meteorology and Climatology , American Elsevier Publishing Co., Inc., New 
York. 

POLLACK, J. B., AND P. BODENHEIMER 1989. Theories of the origin 
and evolution of the giant planets. In Origin and evolution of atmospheres 
(S. K. Atreya, J. B. Pollack and M. S. Matthews, Eds.), 564-602. 

POYNTER, R. L., AND R. K. KAKAR 1975. The microwave frequencies, 
line parameters, and spectral constants for NHj. Astrophys. J. Suppl. 29, 
87-96. 

RATHER, J. D. G., B. L. UlICH, AND P. A. R. ADE 1974. Planetary 
brightness temperature measurements at 1.4 mm wavelength. Icarus 23, 
448-453. 

ROMANI, P. N., I. DEPATER, AND S. K. ATREYA 1989. Neptune’s deep 
atmosphere revealed. Geophys. Res. Lett. 16, 933-936. 

Rudy, D. J., D. O. Muhlmenam, G. L. Berge, B. M. Jakosky, 
P. R. CHRISTENSEN 1987. Mars: VLA observations of the northern 
hemisphere and the northern polar region at wavelengths of 2 and 6 cm. 
Icarus 71, 159-177. 

SILL, G., U. Fink, AND J. R. FERRARO 1980. Absorption coefficients of 
solid NH S from 50-7000 cm -1 . J. Opt. Soc. Am. 70, 724-739. 

SPILKER, T. R. 1990. Laboratory measurements of Microwave absorptivity 
and refractivity spectra of gas mixtures applicable to giant planet atmo- 
spheres. Ph. D. thesis, Stanford University, Stanford. 

STEFFES, P. G., AND J. M. JENKINS 1987. Laboratory measurements of 
the microwave opacity of gaseous ammonia under simulated conditions for 
the Jovian atmosphere. Icarus 72, 35-47. 

SUTTON, E. C., G. A. Blake, C. R. Masson, and T. G. Phillips 
1985. Molecular line survey of Orion A from 215 to 247 GHz. Astrophys. 
J. Suppl. 58, 341-378. 



138 


TERMAN, F. E. 1943. Radio Engineers’ Handbook. McGraw-Hill Book Com- 
pany, Inc., New York. 

TOWNES, C. H., AND A. L. SHAWLOW 1955. Microwave Spectroscopy. 
McGraw-Hill, New York. 

TYLER, G. L., D. N. SWEETNAM, J. D. ANDERSON, S. E. BORUTZKI, 
J. K. Campbell, V. R. Eshleman, D. L. Gresh, E. M. Gur- 
ROLA, D. P. HINSON, N. KAWASHIMA, E. R. KURSINSKI, G. S. LEVY, 
G. F. Lindal, J. R. Lyons, E. A. Marouf, P. A. Rosen, 
R. A. Simpson, AND G. E. Wood 1989. Voyager radio science obser- 
vations of Neptune and Triton. Science 86, 1466-1473. 

ULABY, F. T., R. K. MOORE, A. K. FUNG 1981. Microwave Re- 
mote Sensing Active and Passive. Addison-Wesley Publishing Company, 
Inc., Reading, MA. 

ULICH, B. L. 1974. Absolute brightness temperature measurements at 2.1 
mm wavelength. Icarus 21, 254-261. 

ULICH, B. L. 1981. Millimeter-wavelength continuum calibration sources. 
Astron. J. 86, 1619-1626. 

ULICH, B. L., J. R. COGDELL, AND J. H. DAVIS 1973. Planetary 
brightness temperature measurements at 8.6mm and 3.1 mm wavelengths. 
Icarus 19, 59-82. 

ULICH, B. L., J. H. Davis, P. J. Rhodes, and J. M. Hollis 
1980. Absolute brightness temperature measurements at 3.5 mm wave- 
length. IEEE Trans. Ant. Prop., AP-28, No. 3, 367-377. 

ULICH, B. L., J. R. DICKEL, AND I. DEPATER 1984. Planetary Obser- 
vations at a wavelength of 1.32 mm. Icarus 60, 590-598. 

VALKENBURG, E. P. AND V. E. DERR 1966. A High-Q Fabry-Perot 
interferometer for water vapor measurements in the 100 Gc/s to 300 Gc/s 
frequency range. Proc. IEEE 54, No. 4., 493-498. 

VAN VLECK, J. H., AND V. F. WEISSKOPF 1945. On the shape of 
collision-broadened lines. Rev. Mod. Phys. 17, 433-443. 

WALLACE, L. 1980. The structure of the Uranus atmosphere. Icarus 43, 
231-259. 



139 


WATERS, J. W. 1976 Absorption and emission of microwave radiation by 
atmospheric gases. In Methods of Experimental Physics, M. L. Meeks, ed., 
12 Part B, Radio Astronomy, Academic Press, Section 2.3. 

Werner M. W., G. Neugebauer, J. R. Houck, and M. G. Hauser 
1978. One millimeter brightness temperatures of the planets. Icarus 35, 
289-296. 

WEST R. A., D. F. STROBEL, AND M. G. Tomasko 1986. Clouds, 
aerosols, and photochemistry in the Jovian atmosphere. Icarus 65, 161-217. 

WEIDENSCHILLING, S. J., AND J. S. Lewis 1973. Atmospheric and cloud 
structures of the Jovian planets. Icarus 20, 465-476. 

Willey, D. R., T. M. Goyette, W. L. Ebenstein, D. N. Bittner, 
AND F. C. DeLUCIA 1989. Collisionally cooled spectroscopy: Pressure 
broadening below 5K. J. Chem. Phys. 91, 122-125. 

ZHARKOV, V. N. AND V. P. Trubitsyn 1978. Physics of Planetary 
Interiors (W. B. Hubbard, Ed.), vol. 6 of Astronomy and Astrophysics Series. 
Fuchart Publishing House, Tucson, AZ. 

ZHEVAKIN, S. A., AND A. P. Naumov 1963. Coefficient of absorption of 
electromagnetic waves by water vapor in the range 10 microns-2 cm. Radio - 
phys. Quantum Electron. 6, 674-694. 



